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Crossing Symmetry and the Relativistic Equation 
for two Fermions+ 


By 8. F. Epwarps and P. T. Marruewst 
Department of Mathematical Physics, University of Birmingham 


[Received March 19, 1957] 


ABSTRACT 


A consideration of the symmetries of the scattering matrix of two 
fermions leads to the proposal of a linear scattering equation. This 
equation, unlike the ladder approximation, is shown to have all the 
symmetries and invariance properties of the true ¢ matrix. In the limit 
of one of the fermions becoming indefinitely heavy the relation of this: 
and the ladder approximations to the Dirac equation is discussed and it 
is shown how a covariant equation which reduces to the Dirac equation: 
in the no-recoil can be formulated. 


§ 1. LyTRODUCTION 


THE construction of relativistic two body scattering equations is a well- 
known procedure from the work of Bethe and Salpeter (1951), and 
Schwinger (1951). In this formulation the equation for the two body 
propagator (or. wave function) is a clear generalization of the non- 
relativistic equation, and can also be viewed as a summation of a series 
of Feynman integral, in the simplest case the ladder summation. How- 
ever, the complete solution of the problem satisfies several invariance 
requirements, in particular a relativistic theory must exhibit the correct 
statistics between ingoing and outgoing particles as well as the usual 
symmetries amongst ingoing and outgoing where appropriate. The 
former crossing symmetry is a very powerful restriction in the theory of 
nucleon—meson scattering and must necessarily be satisfied by any theory 
of this reaction (Chew and Low 1956, Feldman and Matthews 1956). In: 
particular, if one wishes to use the relativistic analogue of ordinary 
potential scattering, the Bethe-Salpeter equation, which in the ladder 
approximation violates this symmetry, has to be modified to include it 
before it can be successfully applied. After modification a very simple: 
theory of meson—nucleon scattering results, which agrees with experiment, 
and with discussions based upon unitary and causality which are sufficient 
alone in no-recoil theory to indicate the form the solution will take (Chew 
and Low 1956, Edwards and Matthews 1957 a). This paper discusses 
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whether the crossing symmetry imposes any serious restrictions upon 
processes involving two fermions and whether the theory can be modified 
to admit this symmetry. It has been noted elsewhere (Edwards and 
Matthews 1957, henceforth I) that the crossing symmetry does indeed 
exist for two fermion scattering, but is of vital importance only if one of 
them can annihilate with the other. For this case the modified equation 
for the propagator has to be used before a successful theory can be 
evolved. In the limit of the mass of one of the fermions tending to 
infinity, the exact equation for the other is of course the Dirac equation, 
but this is not easily obtained from the relativistic equations. This 
difficulty is discussed and it is shown how the modified equation improves 
the situation. 


§ 2. CrossInG SYMMETRY FOR Two FERMION SCATTERING 
The crossing symmetry for fermions has been discussed in detail 
elsewhere (I), so here we shall simply quote the results. If the process is 
illustrated diagrammatically as in fig. 1 with p,, p, incoming momenta and 
945 Jz Outgoing momenta, considering the particles as distinguishable for 
the moment, then if one introduces 


N=Pivn 

Te=PotWe “18s Sel ape eR 
fe Pi == Feet 
the ¢ matrix satisfies 

[e(r,, 9, ls; yD, yO"); =t (Sr, Yo; 3; yh a peed) J :. (2) 
=AMO(ry, 19, Po YD, yO"), (8) 
The indices Tr(1), Tr(2) mean transpose of all matrices in the 1 and 2 
spaces respectively. This algebraic operation is equivalent to swinging 


over the | lines, then reversing the arrow on the line (fig. 2), and similarly 
for the 2 line. 


Fig. 1 Fig. 2 


PA UU 1! pi 


pa : p2 1 


These remarks are true only for even interactions, i.e. those which do 
not change sign under charge conjugation, e.g. the pseudoscalar meson 
interaction y;. For interactions involving a change of sign under charge 
conjugation the ¢ matrix has to be considered a function of the coupling 
constant which changes sign upon crossing. The diagrams illustrated 
are for particle—particle scattering but of course the ¢ matrix also includes 


and the Relativistic Equation for Two Fermions 833 


particle anti-particle scattering which diagrammatically amounts to a 
reversal of arrows on the anti-particle line. 

In the centre of mass system, the meson—nucleon crossing symmetry 
operation transforms the energy w and cosine of the angle between 
initial and final relative momentum z(=h . k1/k?, k2 =k!) into w, and 2, 
where 

wy V/ (m+ wP?—p?)+3(w?—p?)(142,) 


=—w/(m+o%—p)—Ho®—p(If2) . 2 2. (4) 
(w®@—p2)\(1—z)=(wy@—p2(I—z,). 8) 

For m>w these become 
w= —w+(1+2)(w®—p2)(m-+ 20)! 00m)... (6) 
24=2—2m(1—z*)(m—2wz)-1+ 0(m-3) ee ae a 


and so correspond approximately to reflection in energy. It is shown in 
I that for fermion—fermion scattering the corresponding relation is 


wP=—(w?—p?)(1+2)/2. pote en ee |S) 


Whereas in the meson—nucleon case the ¢ matrix at w is approximately 
related to its value at w, which by simple manipulations can, for example 
in effective range theory, be related back to the value at w, the form (8) 
offers little information on the structure of the corresponding ¢ matrix, 
since there is no simple relation between w, and w or indeed any physical 
energy. However, when say nucleon—anti-nucleon scattering is con- 
sidered, it can be shown (I) that the scattering contains a part t,, which 
satisfies (8) and can be expressed in terms of nucleon-nucleon scattering 
(except for states for which the latter is excluded by statistics), and also a 
part ¢, which satisfies the crossing relation 


W'==W4", 24==—2. eae ee on wan Ce) 


This is a restriction as powerful as that of meson—nucleon scattering 
and any theory must conform to it. Even in the nucleon-nucleon case 
the observance of the crossing symmetry brings the theory much nearer 
to a form which allows non-relativistic and no-recoil limits to be taken 
correctly. In § 3 the theory will be set up in invariant form. 


§ 3. THE INVARIANT THEORY 


In this section is derived the simplest relativistic generalization of the 
non-relativistic theory which satisfies all the invariance properties dis- 
cussed above. This will neglect self energy effects, i.e. virtual emissions 
and absorption in which the same fermion emits as absorbs, though the 
formalism can be extended without difficulty to these processes if required. 
A rigorous derivation of the relativistic equation is given in the Appendix 
using Schwinger’s functional differential formalism, whilst: in what 
follows we shall put it forward as a plausible consequence of our invariance 
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requirements. Consider for example two nucleons (taken not identical) ; 
the ladder approximation for their propagator is 


(94-25 y+. 12)Gaa(1, 2 33, 4)=38(1—2)8(3—4) . . . (10) 
eye tad 11 
-where iS =F 0x0 + M4): 6 : 5 3 5 P ( ) 


Similarly S,-1! and J,, is the interaction which in ladder approximation 
of pseudoscalar meson theory is g?y;5')y,5'?)4(13)6(2—4)6(3—4). The ¢ 
matrix is obtained from G,, by the substitution 


G1 5=S,8.+8,St5185, . . . . . ° (12) 
integrations being suppressed, and satisfies 
dele 2134) Teo 34) 
+ fTsa(l, 5 ; 36)S,(5, 7)S_(6, 8)t,(7, 2; 8, 4) 
<d¥ 1dted*6d46. =. 2%, 3 3) a. Se 


This equation is shown symbolically in fig. 3. It contains the series of 
ladder diagrams when the first approximation for J is taken, and omits 
crossed diagrams, i.e. those in which meson lines cross. In the Bethe- 
Salpeter approach these diagrams are included by altering J,,, and in 
principle J,, can be written down to make the equation exact, the formal 
expression for example in terms of functional derivatives having been 


Fig. 3 


a 


given by Schwinger (1951). Here it is proposed to symmetrize the kernel 
of (13) by adding the crossed term derived. from it (fig. 4). To have the 
crossing also on the lower fermion line (i.e. time reversal) the whole thing 
should be symmetrized, leading finally to 
Peres =— 71, 34)-+4f | [f d45d*6d!7d48 

X [Tia(1, 5 ; 8, 6)9,(5, 7)S,(6, 8)£,.(7, 2 ; 8, 4) 

si ya(2;'0 5 3, L)Ss(De S108 )ias( Loleeoe se) 


and the Relativistic Equation for Two Fermions 835 


which is given diagrammatically in fig. 5. J in this equation now satisfies 
all the symmetry properties of ¢ in all approximations as it clearly does 
in the first. For the complete analysis both (13) and (14) can be exact 
(with different J function of course), but (14) satisfies all the symmetry 
and invariance properties for ¢ in its approximate form. The original 
form (10) can be discussed as a homogeneous equation, i.e. a wave equation, 
but this is not possible in the symmetrical form. The reason is that we are 
emphasizing the symmetry relations between initial and final states, 
which is not an easy concept in terms of wave functions. 


I 
k 


Using Born approximation for J the new equation fails at the third 
iteration since it gives the diagrams of fig. 5 with a coefficient of } instead 
of unity. This has to be compensated for by a change in J without altering 
the invariance property of J. 

When the two fermions can interchange through the interaction, e.g. 
proton, neutron and charged meson, / will contain the two types of 
diagram of fig. 6. For nucleon—anti-nucleon scattering these terms are 
given by fig. 7. 
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This structure of terms containing a pair annihilation and creation, 
and those in which nucleon and anti-nucleon lines run continuously 
through a graph, typify the two classes of term which the diagram series 
can be split into. Though this is not an easy concept mathematically, the 
diagrams reveal very simply the fact that part of the nucleon—anti-nucleon 
scattering resembles nucleon-nucleon, whereas part is quite different. 


§ 4. THE No-Recoi Liwrr 


If one of the fermions becomes infinitely heavy, the problem is exactly 
determined by the Dirac equation for the light fermion in the potential 
of the heavy. Being exact the Dirac equation must be the sum of ali the 
diagrams and it is easy to see how this comes about. The propagators 
of the heavy fermions become step functions in time and delta functions 
of space 

S(1, 2) > O(t, t,) exp[(im(t,—t,)]6(w7;—w.) . . . (15) 
and an integration over boson lines entering the heavy fermion line, when 
all diagrams are summed, is just an integration of the time coordinate of 
the vertex from — oo to + «, all the step functions adding up to unity. 
For example, in second order we have fig. 8 (parts of diagrams not shown 
being identical). 


These give, when the exponential factors are removed by end conditions, 
the integrals 


I fe dt,dts[ (A, 1)A(B, 2)O(1, 2)-+A(A, 1)A(B, 2)(2, 1)] 


={_ { diylt,A(A1)A(A2) = V(a)V(b) ss. (16) 


where a, 6 are the (space) distances of A, B from the origin where the 
heavy fermion is, and V(a) the static potential 


Via)= [< A(A, Udt,. 


Since all diagrams are required to get the no-recoil limit it is clear that no 
approximation which can be expressed in terms of diagrams can ever 
give the correct limit, though clearly eqn. (14) is very much closer to it 
than eqn. (13) which does not even resemble the symmetry required to 
take the limit. : 

Since no ‘set of diagrams’ can give the required limit it is interesting 
to note that there is a simple addition to (14) (though not a unique 
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addition), which gives an estimate of the missing term, is covariant and 
yields the correct no-recoil limit. Say for example we are considering 
electron—proton scattering given diagrammatically by fig. 9, 
Ey ae i ea iL, 3)5(1—2)8(3—4) 

+4y,My,D(I, 3) {fsa B)S(1, 6)t(5, 6; 3, 4)d*5d26 

-+sym. terms. | 
Then the addition of the term 


1 supe ; ni i : 
a | {7 PQ, 5)S8,(, 6) {74° {16, 2; 3, 4), (asa apm pp 4) 
x d45d*6 
provides a possible extension. (The brackets { } stand here for anti- 


communicator.) This can be thought of as a ‘diagram’ as in fig. 10. 
Indeed in the field theoretic derivation of the Appendix it is represented 


Fig. 9 


4- SY mM. 
“ Cerngs 


by a functional derivative of the ¢ matrix. This term removes an So, 
and replaces it by an S,yS, as required. If this procedure were exact 
this form would be symmetric between the two fermions, and if this 
symmetry is important clearly a term symmetric between the two 
fermions could be employed. The S function in the denominator could 
of course never occur in the more usual kind of expansion. 
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The Derivation of the Equation 


This most convenient technique for the calculation of equations of this 
type is the functional differential, introduced by Schwinger. Since we 
deal only with nucleon-nucleon scattering and shall not consider radiative 
corrections it is convenient to eliminate right away the meson variables, 
whereupon the effective equations of motion become 


(20,4 m—o2 | Acenbln. yobloMy) Houle) (A) 


where 7 is the spinor source. Taking expectation values and differenti- 
ating with respect to 7 
8’ (, «")=8 (ba) )/8y(2") (A2) 


5(K)/8q=i[ (Kb), (K)Cb)] (A3) 
so that 


(4 a,tn—g | Ae, BY, 2), 6A, 2)))8'=Bee—2") (A4) 


Dy, )=8'(y, 2)(8/8 UO) ))— C2): (A5) 
Suppressing integration variables this becomes 
S'=S+9784|Foy.9 718 (A6) 


Since we wish to treat our equation symmetrically we use the symmetric 
form 


S'=S+i9SAJ, yD)S' +4978’ A[Z, ysDIS (A7) 
and in the latter term Y operates from right to left. Now to obtain the 
# matrix, the substitution 


S(wy)tly, © 5 2 )S(z, &")= (8/3 Ceb(£) ))(8/3 b(n) ))S' (ee, a”) (A8) 


“14/5. )SCP))S'S. (A9) 
Applying this to (A6) and multiplying out and performing the function 
derivatives, one obtains 


t=g? A+ t0/( ASSt)+¢? (AS {52/8 (ys 5 ab) Ht) + (A10) 
where c¥ stands for the symmetric combination of (A6), and R for radiative 


corrections. The remaining term, the derivative of t, is dropped in (14) 
and it is essentially this term which is considered in § 4. 


1.e. 
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ABSTRACT 


Some relations between different partial wave amplitudes due to 
‘crossing symmetry are derived for meson—nucleon scattering and some 
similar symmetry properties for the nucleon-nucleon interaction. 


‘ §1. LyrrRopuctTIon 


‘THe Chew—Low (1956) theory of meson—nucleon scattering is based on 
applying the restrictions of unitarity and ‘ crossing ’ symmetry to a fixed- 
source interaction. Their equations have since been rederived by Oehme 
(1956) as the low-energy limit of dispersion relations, which also include 
the limitations imposed by causality. In their exact form the dispersion 
relations are identities, which can only be used as a check on the correct 
interpretation of the physical data. (By interpretation we mean the 
deduction of the scattering amplitude from the observed cross section.) 
It is only after approximation, introducing new physical information, 
that they become genuine equations for the scattering amplitudes. 
‘These equations are very powerful, but they always involve an integral 
over the whole energy range. To evaluate these integrals, or even to 
show that they converge, additional assumptions have to be made. 

We wish to point out here, that crossing symmetry alone implies 
certain relations between the different partial wave amplitudes, which 
do not involve any integrations, and must be satisfied exactly by any 
theory of meson—nucleon scattering. In particular a relation between 
the s- and p-waves is established which indicates that the splitting between 
the 31 and the 13 phase shifts is of the same order as the phase shifts 
themselves. 

We also show that a similar crossing symmetry exists for the scattering 
amplitude of two distinguishable fermions linked by an interaction which 
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is even under charge conjugation. The implications of insisting on this 
invariance in an approximate equation for the two nucleons will be 
discussed in a separate paper. 


§ 2. Crossine In Pron—NUCLEON SCATTERING 


Suppose the meson momenta and charge label are , 7 finally and q,, s 
initially, and the nucleon momenta are p, finally and q, initially. Then 
crossing symmetry, which is a direct consequence of relativistic invariance 
and the Bose statistics of the mesons (Feldman and Matthews 1956). 
implies that for the scattering amplitude 


ts(Pr Px V1: Yo) =bs-(—G1, Po —Pr: 2): Seaway ae (2.1) 
On the energy shell 


Prt Pe=G1 + Ge . ° e 5 = « . (2.2) 
and 


pY=ar =e", pP=qP=m. . . . . . (2.3): 
Consequently there are only three independent 4-vectors and two inde- 
pendent scalar products. For discussion of crossing it is convenient to 
use the variables 


r= (pi+4;)/2 
o=(Pots)/2, Sat eee rome (oe 
3=(P1—41)/2=(P2—o)/2, . 

and to take as the independent scalars 


fVo=G, "F,7==0." 0. Pe. Cee 
Under the crossing operation 
Cae ji—@y 0 05 to ee en eG) 


In the centre of mass coordinate system 
ry=e, —(k+k’)/2, 
P=, (k-EK')/2 ob ee) 
Yg—0, (k’+k)/2, 
where 2k and 2k’ are the relative momenta in the final and initial states 


? 


({k|=|k’|), and « and w are the corresponding nucleon and meson energies. 
Then 

a=2ew+k?(1+-2), 

nes zh i a2) 
where 

Briel 

aa ae 


The scattering amplitude may be expressed in terms of these variables ? 
and according to (1) 


t.(@, z)=t,(w4, 24), ° °. e e e . ° (2.9) 
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where, by (6) and (8) 
2w(m? + w*—p?) 1/2 + (w?—p?)(1+2) 
=~ [20 (m+? —p2) 19+ (op? —w\(1-+2y)) 


and 

 (w®—p2)(1-2)=(02—w)(I—)). . . . (2.10) 

In the limit m— oo, these equations are satisfied by 
pO eee ere 2, el ee) (ANT) 
We thus expand in powers of as and obtain the approximate solutions 
Oe le a aN 1 ste) Kgs (2.12) 

Dad 2w(1—z?) es 

shee Ly Pera (m ip (2.13) 


§ 3. THe PartiaL WAVE RELATIONS 
The isotopic spin projection operators are 
Pr, s=3 Ors t3lt Tel} 
Ps, s=3125,5—31 7p To}: 


ESP a where 
If—1 4 
a= 3 ( 2 1 ) Se Vs OP ae (3.1) 


The angular momentum projection operators are 


Hence 


Si= aA [(2 +1) ¥,(z)—t0,(1—2?) Y;'(2)], 
(3.2) 
S,-= I (ZY, (z z)+t0,,(1—z?) Y/'(z)1 
and 
F 2j7-+1 
Tr Sij(2)S* yy (2) dz= ai : 8195375 (3.3) 


where 77 means the trace of the oc matrices. Then, expanding t in partial 
waves and using (2.9), 


tw, z)=P'S,(z)t; (w), 


= Ptatks*(2,)t,*(w). ea ee Pee e) 
Thus 
8;(z)t;'(w)=a"*s+(2,)t;"(w). cththors AN ee 
By (3.3) to extract the j-wave part of this equation one must operate 
with 
2 2 
| Pr|sy5*(@) eh ae ges oss) (8.0) 
giving 


(21-1)? i ea. 
tyeo)= een) Put Clan 
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Since both z, and w, are functions of z and w, the right-hand side of this. 
equation involves a mixture of partial waves. If t(wj) is expanded as a 
Taylor series about the point w;=—w, the integral may be evaluated. 

At energies up to the resonance one may assume that only S and P 
waves are contributing. Let t(w), (j=1, 3), be the two s-wave amplitudes. 
and #,/(w), (t, j=1, 3) be the four p-wave amplitudes. Then for the 
S-wave equation one obtained 


k? ot = =4wH Rar otie 
Of (a — pik} zk 0 Oy EE DO a Ne ee el eee ae ‘ 
t'(—w)=a E (w) aan 37, 1 rece +0O(m )| : ESS 
where TT*=1t,"+ 2t,*, No wee ee a a 
and for the p-wave equation 
Pot areal 
Of Oa —gikg. k Sg NE MRL LE Le D ™ =o ; 
t;'(—w)=a""Gjn, tn (w) wee Rie Fee +O(m )|. :. (3210) 


where a,,, is the matrix (3.1), but now operates on the angular momentum 


variable and 
Gaal Le as ek ore eee eo 


§ 4. Discussion 


Equations (3.8) and (3.10) are our main results. Taking only the first 
term on the right-hand side of (3.10) gives the relation between the various 
p-waves in the static limit, which has been obtained previously by Chew 
and Low (1956). The additional term shows a further dependence upon 
the S-waves. If one subtracts the equation for t,° from that for t,! one 
obtains 

Kee 
t+) )= 9 ae [taw)—t*(w)+E(w)—P(w)], . (4.1) 
where ¢(-++) is the even part of #(w). If the S-waves are neglected, this. 
is clearly consistent with the well-known degeneracy 


Pa te . . . . : - - . (4.2) 
According to Drell et al. (1956) the experimental points may be 
approximately fitted by 


t1(w)= —0.04+ 0.140, 
4.3) 
£3(w)=—0.04—0.07u, We 
in units of uw. Taking these values 
ke 
pel Sy 3 FNS Re ‘ 
ts'(-+)—t,3(-+) Smelt + ss. (4d) 


which indicates a splitting of about 1° at ke. This is of the same order 
of magnitude as the phase shifts themselves, and shows that these two. 
phase shifts must be treated separately and in conjunction with the s-waves 
in any adequate theory. 

Equation (3.8) shows how the s-waves are coupled through crossing 
symmetry to the derivative of the p-wave. y 
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This investigation was originally undertaken in the hope of relating the 
energy dependence of the s-waves to the resonance in the (3.3) state. 
We have been unable to find any indication of this. 


§ 5. FeRMIOoN—FERMION CRossinG SYMMETRY 


Consider the interaction of two distinguishable fermions, through an 
interaction which is even under charge conjugation. Then the propa- 
gation corresponding to scattering is 


G(D1, Pos d1 W2)= | exp [—2(p % + Pol o—G1Y1—FeY 2) | | (5.1) 
CD (b(@y)0(a2)1(Y1)ba(Yo))) Try... dy, 
=| exp [—0(p 4% 4+ Pol o—G1Y1— 122) | 


(—TbaYa ales (es balya)) dey... Ty, (5.2) 

Now the charge conjugate field is given by 
vi=—Op, =O, . . . . 2. . (5.3) 
C'=—C, C ly, C=—y,' oe) 


Thus, substituting for %,, in terms of 4,1 and interchanging the dummy 
variables x1, 7, 


Gyp(P1> Po I> 2) = | exp [—1(G4@ 1+ PX 2+ PiY1—Wey 2) | 
Cup T (1 (%1); ho(2o)h4' (Ys)bo(Y2)) C-4g d*x,...d*ys 
=[C771G*)(—9y, Po. —Py Jo)C 1)? Pe peer RS): 
where ¢(1) implies the spinor transpose with respect to particle 1 and G 
denotes that particle 1 has been replaced by its anti-particle. 


However for even interactions the anti-particle interactions are the 
same as the particle interactions. Thus, using (4.4) we have 


(py, Pa Tv V2: yh), io) 
=G'O)(—qy, Pa» —P1r; 2 —y(1)’, y(2)). age! ons (5.6) 
The effect of the double tranpose is to reverse the order of all the (1) 


factors. 
In terms of the variables defined in § 2, this reads 


Gry, V9, "35 y(1), y(2)) 
eC = ea Tene Vlg) eae ct oe (8.7) 


Insisting on this invariance in approximate equations brings the theory 
much nearer to a form which allows non-relativistic and no-recoil limits 
to be taken correctly. This will be discussed in detail in a separate paper. 

After this manuscript was completed the author’s attention was drawn 
to some similar considerations by F. E. Low (1956), the results of which 
were reported at the International Conference in Pisa. 
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SUMMARY 


In this paper the problem of the thermal instability of an incompressible 
rotating fluid sphere heated within is considered. The equations govern- 
ing the state of marginal instability are derived for the case when the 
motions and the associated perturbations have symmetry about the 
axis of rotation. The underlying characteristic value problem is solved 
for a slightly modified set of boundary conditions. Nevertheless, the 
solution obtained suffices to determine the general nature of the depen- 
dence of the lowest Rayleigh number for the onset of instability on the 
Taylor number 7'(=42?R4/p?). 


§ 1. INTRODUCTION 


In recent years the onset of thermal instability in horizontal layers of fluid 
heated from below has been investigated both theoretically and experi- 
mentally under a wide variety of conditions (Chandrasekhar 1952 a, 
1953 a, 1954a, b, 1956a, Chandrasekhar and Elbert 1955, Fultz 
and Nakagawa 1955, Nakagawa 1957). These studies have disclosed a 
number of novel features: the onset of thermal instability as oscillations 
of increasing amplitude (i.e. as overstability) when a metallic liquid such 
as mercury is subject to Coriolis acceleration, and the conflicting nature 
of the effects of rotation and magnetic field when they act simultaneously. 
While these studies have been useful in thus drawing attention to the 
different types of phenomena one may encounter when rotation and 
magnetic field are present, it is important for applications to problems 
such as convection in the earth’s core that the relevant investigations 
are extended to a spherical geometry. The simplest problem in this latter 
category, that of the thermal instability of a fluid sphere heated within, 
has indeed been investigated (e.g. Wasiutynski 1946, Jeffreys and Bland 
1951, Chandrasekhar 1952 b, 1953 b, Backus 1955); but the extension to 
the case of a rotating sphere or when a magnetic field is present raises 
mathematical difficulties of a different order from those encountered in the 


7; Communicated by the Author. 
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solution of the corresponding plane problems. The reasons for the 
increased difficulty are principally two: first, while the physical circum- 
stances, at best, allow only symmetry about an axis, the boundary 
conditions are specified on a sphere; and second, it does not appear that 
the problems can be solved by variational techniques such as those 
which enabled the solution of the plane problems. 

The present paper is devoted to the problem of the thermal instability of 
a rotating fluid sphere with-internal heat sources. The manner in which 
the solution is obtained in one very special case may be suggestive of 
methods one might adopt for the solution of more difficult problems. 


§ 2. THe EQuaTIONS OF THE PROBLEM 


We shall consider a homogeneous sphere of radius #& rotating with an 
angular velocity @ about the z-axis and with a uniform distribution of 
heat sources such that in the absence of conduction the temperature at 
each point will rise at the rate «. In treating this problem we shall ignore 
the rotational flattening of the sphere: this is justified in our present 
context, since the principal effect we are seeking is that of the Coriolis- 
acceleration term, 2ux< 2, in the equation of motion; for this purpose 
the small departure of the equilibrium configuration from a sphere is not 
essential. 

The temperature distribution in the equilibrium state is given by 
(Chandrasekhar 1952 b, eqn. (2): this paper will be referred to hereafter 
as I) 

T=81R—r?), B=e/6x, < ot ie oe et 


where « is the coefficient of thermometric conductivity. Under these 
circumstances, the equation governing small departures from the initial 
stationary state are (cf. I, eqns. (11), (13) and (14)): 


ay 
5, kK V70+2Bu .r Sh a i rT 


and 
du 
at 
where » denotes the coefficient of kinematic viscosity, p the density, 
u the velocity, 6 and dp are the perturbations in the temperature and 
pressure respectively, dw=dp/p and I, is a unit vector along the axis of 
rotation, and 


=— grad da+yér—y curl? u+-2.Qux I, a) ns) 


4 
VY — 3 7Gup. A A , . 5, ° - . (4) 


In eqn (4), G is the constant of gravitation and « is the coefficient of 
volume expansion. In addition we have the condition 

div u=0, + cieietet ric pace eC 
on the velocity field. 
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§ 3. THE Equations GOVERNING MARGINAL STABILITY IN CASE THE 
VeLociry Fretp HAS AxtaAL SYMMETRY 


In this paper we shall restrict our consideration to the case when the 
instability sets in as a stationary pattern of convection. In this case the 
equations governing marginal stability are : 


BY Ca 20U IS eo 4 oo 8) 
and 
grad da=yér—v curl? u+2Qu~ I,. he Se) 


In the further discussion of these equations we shall assume that the 
velocity field as well as the distribution of 6 have symmetry about the 
Z-axis. 

If the velocity field is divergence free and has axial symmetry, we can 
express it as a superposition of a poloidal and a toroidal field in terms of 
two sealars ; thus (cf. Liist and Schliiter 1954, Chandrasekhar 1956 b, c) : 


Ussler ie curl (Xr) in ee et) 


where J’ and U are two scalar functions which are azimuth independent. 
The physical interpretation of these two scalars is the following: U 
defines motions which are entirely in the meridional planes and V defines 
motions which are entirely rotational. 

The particular advantage of the representation (8) in the present 
connection is that we can write down without any difficulty the result of 
the operation, any number of times, of curlonu. Thus (cf. Chandrasekhar 
1956 b, eqn. (18)) 

eurl u=—I,xrd,U-+ curl (I,xrV), CO a ED, 


where 4, is the Laplacian operator for axisymmetric functions in‘ five 
dimensional Euclidean space; by repeated applications of this formula 


we obtain 
curl? u= —I,x r4,V —curl (I,xr4,U), ee CLO 


core lord, —curl (IlexrdeVy," 2 2 (14) 


etc. 
Returning to eqns. (6) and (7), we first eliminate 6@ by taking the 


curl of eqn. (7). We thus obtain 


y grad @Xr—v curl? u+22 curl (uxi,)=0. . . . (12) 

It can be readily verified that 
10xr=l fe ag ha (13) 
grad @xr=1,xXr{ a —— a ), 3) 


where w[=1/(x?+y?)] is the distance of a point from the axis of sym- 
metry; also, 


oU 
curl (ax Lat xr Se teu (Lx 5). CAL 


3M2 
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$45 
Making use of eqns. (11), (13) and (14), we can rewrite, eqn. (12) in the 


form 


06 =z OO 


) rfl. xrd,tU—eul (1, x47) 


70) eee OU em 
Sager az +eourl(I,xr Oe =, 
The poloidal and the toroidal fields represented in this equation must 


separately vanish. Thus, we must have 
é ot x)? (16) 


(15) 


2Q20Vy 
DY fee be ond 
ne v Oz v\0z2 wdw 


220U 


— =0. 


and 
AsV+ y Oz 


Now applying the operator 4; to eqn. (16), we have 
0 0 
( ‘ ) ea i aa eee) 


220 y 
3 pe — i ce Cee 
ks cosy oe acammviene Vase see 

since the vperations 4, and d/dz permute. Eliminating V from this 
last equation by making use of eqn. (17), we obtain 
40°0?U 240 
3 at ee — ee 

Aaa y= dz" L As 5 ==) eames h), 


7] 
A ei xn erie (Pe a 
dz wo Ow) 0z w0m 
where 4, is the usual three dimensional Laplacian for axisymmetric 


It may be readily verified that 
z 0 Gr geno 
}=( ) Are me 
By making use of eqn. (6), we may therefore rewrite the 


functions. 
right-hand side of eqn. (19) in the manner : 
=) 
Ur. 


On the other hand, according to eqn. (8) 
Ges) 
ura—(5— 22) ov (22) 
(23) 


Equation (19) thus reduces to 
4022 Q2U c3 2By 


02 wou 


Oy eran\2 
AsO eae = el ig) miu. 
In spherical polar coordinates (= cos #, @) eqn. (23) takes the form 
4Q2/ 9  1—p?A\2_  2By a 
a. ————— = ey gees Be. ‘ 
ae Ot y (Hs, r x) L KV Beet hee - (24) 


Ls? 0? = Bind 
ee... (28) 


re uz 72 Om 


where, now 
+ 302 eee 
> Orr Or 
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Similarly, eqns. (16) and (17) take the forms 
He ese Q0. 1—p 
A30— =" (, nts) ya vie ey azo 


O Or r Op vr op’ 3 
and 
eee tO. lug 
4,V+ —{ p+ —-)U= 27) 
ve = (uxt : 5) U=0 eat ia arn 
Also, it may be noted here that 


a (1—p2)1/2 Q 
eee lee oh ee By ft cotile 

leh eT See ove Catal Re athe Ral 
where I,, Ij and I, are unit vectors along the ares dr, rd3 and r sin #d@ in 
the three principal directions. 


Pp? (28) 


3.1. Boundary Conditions 
Solutions of eqns. (24), (26) and (27) must be sought which satisfy 
certain boundary conditions. These depend on the nature of the bounding 
surface at rR; but in all cases 6 and the component of u in the direction 
I, must vanish. According to eqns. (26) and (28) these latter conditions. 
require : 
Se ee pees] 0 al 0 
V==0 and fe U—(u5 + rei Os 
If the bounding surface ;=F is rigid and no slip occurs here, then the 
components ws and u, of u along Igy and |, must also vanish. From 
eqn. (28), it follows that the vanishing of these additional components of 
u at 7=F require 
Oly ty ; at 
Or =—=(Qand /)=0 on r= (rigid surface). . . . . (30): 
(Actually, the vanishing of ws requires only that ¢(r?U)/dr=0; since, 
however U=0 at the boundary, the condition on U stated follows.) 
On the other hand, if the bounding surface should be free, then the 
vanishing of the tangential viscous stresses on this surface requires 


(I, eqn. (21)) 
0 (us 0 (Up 
se el ea a Rapti 


According to eqn. (28) these conditions are equivalent to 


) V=00onr=—Rk. (29): 


Of ld CU e200 20 ] 
Set. == ("7 ) ee a Se Oe 
Or| r2 or or? r or r2 
and ee (Oz) 
= =0 On fai. | 


Since U=0 on r=R, the first of these conditions can also be expressed 
as 


2 


Srv)=0, re seas) 
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Collecting the foregoing conditions, we may write 


J a _ is eS 
== 0. Eh. G 5 V=0 and Ze U) 3 er 3 ( ) 
on a rigid hounding surface, and 
Ge Vs 201—waV __ . 
U=0; 5a U)=0, am (ange 0 ee reer =0, . (35) 


on a free bounding surface. Finally, we may note that if we measure 
distances in units of the radius of the configuration, we may rewrite 
eqns. (24) and (27) in the forms: 


02U Gh 
Bip T= =O (Se | ne 2 
ASU+T 55 =C 5 [0-H] (36) 
and 
KO) - NO! 
ys nly (Pe a5 
As : R Ae? (37) 
where 
2 D 
pleietys and Cee ee er) Kees 
Vv KV 


re the Taylor and the Rayleigh numbers, respectively. 


§4, THE SOLUTION OF THE CHARACTERISTIC VALUE PROBLEM 
FOR A SPECIAL CASE 


It does not appear that the characteristic value problem presented by 
eqns. (24) to (27) together with the boundary conditions (34) or (35) can 
be solved in any simple way. In particular, it does not appear that the 
problem can be solved by a variational method. However, since our 
principal interest is in the asymptotic dependence of C on the Taylor 
number 7’, the following method of obtaining this dependence suggests 
itself. 

From our experience with the plane problems we may feel assured that 
as a general rule the nature of the dependence of the critical Rayleigh 
number on auxiliary parameters such as 7' (to allow for the effects of 
Coriolis acceleration) or Q (to allow for the effects of a magnetic field) 
is largely independent of the boundary conditions with respect to which 
the underlying characteristic value problem may have been solved ; indeed, 
in all cases which have been examined the asymptotic dependences on 
parameters such as 7’ and @ are the same for solutions satisfying very 
different boundary conditions. Consequently, it would appear that we 
may seek solutions for the spherical problems also with the boundary 
conditions slightly modified so as to allow a simpler solution of the 
characteristic value problem; this will be permissible, if our particular 
concern is to discover the nature of the dependence of C on the Taylor 
number, for example. 
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Thus, if we should consider the problem without rotation (treated in I) 
we are required to seek solutions of the equation (in our present notation) 
2 

AgU=C SU, Oks tetas poe A 
together with the boundary conditions 

U=0) 4200" fr=1) 
and 

either oe ==0 or 57a (TU)=0 pa l) ia alen to (AG) 

depending on whether the bounding sphere is rigid or free. Solutions 
appropriate for these boundary conditions have been obtained by vari- 
ational methods in I. However, for the slightly modified boundary 
conditions 

ou. 400 

ore | + Or 
the problem admits of an explicit solution. For (Chandrasekhar 1956 d, 
eqn. (36)) 


=(Q.cand?) 2°00" 2. & (4)) 


ee a ee 


where «, , is the jth zero of the Bessel function J,,,3/.(7) and C,,3/2() 
is the Gegenbauer polynomial, satisfies the boundary conditions (41) 
and represents a solution of eqn. (39) provided 

af; ,=(n+1)(n4+2)C. eI Co een cS) 
If we compare this formula for C with that given in I (eqn. (81); note 
that the order of the spherical harmonic / in this equation is related to the 
order of the Gegenbauer polynomial, n, by /=n-+1) we observe that the 
present solution predicts the same general dependence of C on the order 
m as that obtained in I for solutions satisfying the correct boundary 
conditions; moreover, the asymptotic dependence of C on n for n— oo 
is the same for all three boundary conditions (namely, 


@U/ar=0, d%X(r2U)/ar2=0 and 9%(r4U)/Ar2=0). 


For the reasons stated, we shall seek a solution of the equation 


aU o* 
sU/ Say Pe a Ve . . . . 4 
A3U+T a2 Cra u?)U, (44) 
which satisfies the boundary conditions 
2 
t=), ag ob eo =(0 and @=0 on r=l. “ee yee C259) 
Orr oF 


We can satisfy the boundary conditions on U by expanding it in 
terms of the various modes 7',, ; (which it may be noted here form a 
complete set of orthogonal functions); thus, 

bay AG 


ni 


(46) 


Ope 
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When the characteristic value C has been determined by this substitution 
for U, we shall also have satisfied the boundary condition on @: for, 
according to-eqns. (24), (26) and (27) 
100\ 2B a ae 
NS = (Li) Ure hoe eee 
when U has the form (46), oe equation becomes 
100% 
=——$ YA, (n+1)(m+2)T ee ee 
Ay(; Ope ih) a K > 
and the required ale of this equation is 
_ 2p 
Ee aos oA 3 il 9 
; oa aor z apes Mn+ 2)T q, 
2p 4 J n+3/al%s, nt) a 


Ss Pager ey errs (Ceeramcaartts Ny 9 
ny) Sin) Te 
or 
2B eA d 
fer Ey Se 1c nsslare a) TfL uO (4) . (50) 
nN; j Jiese 


and this clearly vanishes on 7=1. 

Returning to the solution of eqn. (44) together with the boundary 
conditions (43), we first observe that the results of applying the operators 
4,3 and 0?(1—p?)/du to U given by eqn. (46) are: 


ASU== Fab 08, 2 or 

and we 
5 oe Aaa )U] =—DA,, ,(n+1)(n+2)7,, ; ee) 

nj 


The eftect of the operator 0?/dz% on 7, ; is not so simple. By using the 
various recurrence relations satisfied ay the Bessel functions and the 
Gegenbauer polynomials, we can show that 


a 3 
merntit ase & 5 Pee JGR. AC Oe (FH) : * S (53) 
and 
OT 13/2 | 
ae =Ly, (CP? 4 (tu) +My, (r)C??,(u)+N,, ,(r)C82,_ (1), (54) 
where 


4 aoe’ Nae d [SF najal%, al) jie 3/2 at) = 
Pn. 2n+3 E \ p3/2 ) e | he oe 


; eG m+ 2 d Ds 2%; Aa) ‘ J, (a; r) 
ia n+3 | a ( 378 ) t(n- 3) Sear, - (56) 


a= ae 5) Es j (apd jae n el (ce, , ntl 
PB yea in-+2)) ner Ie (2n4-3)(2n-+5) see soe OL) 


Mik eee (Pn 6m Vas,» Tnssials, of) a 
ag (2n-+1)(2n-++5) y3/2 P) ¥ - . ° e ° . (58) 
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— _ (m9 +2)0 sn I nssyal%s, nf) 
r (2n+-1)(2n+-3) pole 


The functions f,, ;, g,,;, L,,; and N,, , have recently been tabulated 
by Donna Elbert (1957). 
Substituting from eqns. (51), (52) and (54) in eqn. (44) we obtain 


N 


n,s=+(n+1) (59) 


Y An (2%, tint lyin +20} esaialts, oF) care, (,) 
nj 


eae ah k m, Ee ee nC ake aig) = Oe (60) 
m, k 


Equating the terms in the Gegenbauer polynomials (?/?2,, we obtain 


(cf. eqn. (58)) 
| ant bn+ I 
DAs, 1-2, n+l \(n-+-2)C— (n+5)(Qntl)* * DI nsanlty, me 


+P R(4,- 2, pL n—2 eos WV n+2, x)? psla—(), : 2 < (61) 


Multiplying this Tat equation by rJ,,.3/2(%; ,7) and Ses from 
0 to 1 and making use of the orthogonality relation, 


1 
i TI 32%, nT) n3/o(% k, nt) dr=3[J’ n+3/2(%5, Al 35, pce G2) 


we obtain 
: 2n?+6n+1 
[J n+3/2(%;, nA wes, oe ea n (Q2n+5)(2n-+1) 


ee Ase 2: k|n;j)+Anse, (m+2;h|[n;7)}=0 . (63) 


where 


T+(n+1)(n+2)0 5 


1 ‘| 
(n—2;k|n =| P52T xd nag/2l%, nt) a, 
0 


and (64) 


1 
(n+2;k|n; =|, TPN nie, rt nssial%s, nt?) a, 


denote the matrix elements which have to be evaluated numerically. 
Equation (63) leads to the following determinantal equation for the 
characteristic root C’: 
T{(n—2;k|n39)+(n+2;k | 59} 
—$3J nisal%s, aha ree Re Fee 


2n?+6n+1 
x 1%, ates w (2n-+-5)(2n-+1) 


—(n+ n+2)¢! |-« (65) 


From eqn. (65) it follows that 


i on2-+6n-+1 fh 
(n+1)(n+2)C —> 08, fat (Qn-+5)(2n+1) a5, | (T' + 0). (66) 
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A method of solving the infinite order characteristic equation which 
eqn. (65) provides for C would be to set the determinant formed by the 
first few rows and columns of the matrix on the left-hand side of eqn. (65) 
equal to zero and let the order of the determinant take increasingly 
larger values. 


Table 2. The Characteristic Numbers C for the Onset of the Lowest 
Mode of Instability for Various Values of the Taylor Number 

ihe Ist Approx. 2nd Approx. 4th Approx. 6th Approx. 

0 4-116 x 103 4-116 x 108 4-116 x 10° — 

LOOO 6-135 x 108 6-093 x 108 6-056 x 108 — 

2000 8-154 x 108 7-988 x 108 7:846 x 10? — 

4000 1-219 x 104 1-155 x 104 1-105 x 104 — 
6000 1-623 x 104 1-490 x 104 1-389 x 104 1-388 x 104 
8000 2-027 x 104 1-811 « 104 1-650 x 104 1-648 « 104 
10000 = 2-124 x 104 1-894 x 104 1-892 x 104 
20000 — — 2-970 x 104 2-962 x 104 
50000 — — 5-583 x 104 5-557 x 104 
70000 — — ToLST Xe LOF 7-098 x 10 
100000 == — 337156108 9-313 x 104 
200000 — — 1-655 x 10° 1-642 « 10° 
— — 7-267 x 10° Tel 06s LO 


} 1000000 


From eqn. (65) it is clear that the solutions separate into two non- 
combining groups: those which include the even and those which include 
the odd orders of the Gegenbauer polynomials. By restricting oneself 
to determinants of order less than or equal to six by including the terms 
of orders n=0 and 2 in the Gegenbauer polynomials and the first three 
zeros of the Bessel functions J5,. (w) and J,/, (w) it appears that one can 
determine the lowest characteristic root for C with sufficient precision 
for 7’<10°. The secular matrix including the elements for n=0, 2 
and j=1, 2 and 3 is given in table 1. In table 2, the roots determined 
in the manner described are listed ; the corresponding coefficients 4,,; in 
the expansion for U are given in table 3. An examination of this table 
indicates that the values of C have been determined to an accuracy of 
the order of 1%. The results are further illustrated in the figure. 

It is apparent from the figure and table 2 that rotation can be a decisive 
factor in inhibiting thermal convection in a fluid sphere. Indeed, the 
numerical evidence points to the proportionality 


(=—eonstant: Man lias oy 3 wt} (87) 
If this suggested asymptotic relation should be confirmed then (cf. eqns. 


(37) and (38)) 
K [22 ig 
Bo= constant — (5) CEST) iia se aie eels) 
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It is important to remember in this connection that in the presence 


of rotation convection currents cannot be entirely meridional. 


Thus, 


Table 3. The Coefficients A, , in the Expansion U =) An Leas for 
the Solutions Obtained in the Different Approximations (Ay ,=1) 


vi Approx. NG Ao3 Agi Ass As 
1000 second ¥s +0-0288 = 
fourth —0-00010 +0-0286 | +-0-0074 
2000 second es +0-0574 ae 
fourth —0-00039 +0-0561 | +0-0148 
4000 second ay +0-1104 — 
fourth —0-00118 +0-1015 | +0-0274 
second +0-1535 
6000 fourth —0-00175 +0-1322 | +0-0392 
sixth —0-00291 | +0-00032 | +0-1320 | +0-0391 | —0-0019 
8000 second .% +0-1869 
fourth —0-00180 +0-1512 | +.0-0498 rs 
second +0:-2127 
10000 fourth —(0-00133 +0-1627 | +0-0593 
sixth —0-00424 | +0-00042 | +0-1626 | +0-0593 | —0-0032 
second +0-2808 
20000 fourth +-0-00552 +-0-1759 | +0-0957 
sixth —0-00216 | +0-00020 | +0-1758 | +0-0954 | —0-0056 
50000 fourth +0-0279 +0-1605 | +0-1483 
sixth +0-0088 | —0:0041 | +0-1608 | +-0-1474 | —0-010 
70000 fourth +0-0373 +0-1534 | +0-1648 
sixth +0-0126 | —0-0065 | +0-1539 | +0-1636 | —0-0116 
100000 fourth -+0-0463 +0:1469 | +0-1794 
sixth +0-0161 | —0-0094 | +0-1476 | +-0-1779 | —0-0135 
200000 fourth +0-0597 +0-1381 | +0-1995 
sixth +0-0181 | —0-0151 | +0-1392 | +.0-1972 | —0-0172 
1000000 fourth +0-0732 +0-1301 | +0-2183 
sixth +0-011 —0-0237 | +0-1319 | +0-2148 | —0-024 


ne eM | 


according to eqn. (17) rotational motions must be present with an aplitude 
and a character which are determined directly by the meridional motions. 
Finally, it is worth noticing that in the expansion of U in the modes 


Thermal Instability of a Rotating Fluid Sphere 857 


ge (Cl. ean, (42)) the principal terms are the (O71); (2,1) and. (2, 9) 
modes (see table 3). 


6 


logigC 


3 


<5 
4 5 log,oT 6 


Illustrating the dependence of the Rayleigh number C for the onset of the 
lowest mode of instability on the Taylor number 7. 


§ 5. ConcLupING REMARKS 


This paper makes hardly more than a beginning in the study of the 
many problems of hydrodynamic and hydromagnetic stability in rotating 
fluid spheres. Thus, even in the special problem which has been considered, 
the assumption of a stationary marginal state should be given up: for, 
it is known (Chandrasekhar and Elbert 1955, Fultz and N akagawa 1955) 
that in fluids with low Prandtl numbers (as most metallic liquids are) 
instability is more likely to set in as overstability than as convection. 
The discussion of the overstable case was difficult enough in the plane 
problem ; but its importance would justify the efforts which would be 
needed to solve the same problem in a sphere. 


ACKNOWLEDGMENTS 


In conclusion I should like to record my grateful appreciation of the 
assistance I have received from Miss Donna Elbert: she carried out the 
extremely laborious calculations which were needed for this investigation. 

The research reported in this paper has in part been supported by the 
Geophysics Research Directorate of the Air Force Cambridge Research 
Center, Air Research and Development Command, under Contract 
AF19(604)—2046 with the University of Chicago. 


858 On the Thermal Instability of a Rotating Fluid Sphere 


REFERENCES 


Backus, G. E., 1955, Phil. Mag., 46, 1310. 

CHANDRASEKHAR, S., 1952a, Phil. Mag., 48, 501; 1952 b, Ibid., 48, 1317 ; 
1953 a, Proc. roy. Soc. A, 217, 306; 1953 b, Phil. Mag., 44, 233, 1129 ; 
1954.a, [bid., 45, 1177; 19546, Proc. roy Soc. A, 225, 173; 19564, 
Ibid., 287, 476; 1956b, Astrophys. J., 124, 232; 1956c, Proc. nat. 
Acad., 42, 1; 1956d, Astrophys. J., 124, 244. 

OHANDRASEKHAR, oe and ExvBert, D. D., 1955, Proc. roy. Soc. A, 231, 198. 

Exvsert, D. D., 1957, Astrophys. J. Suppl., Series 3 (in press). 

Fuutz, De and NAKAGAWA, Y., 1955, Proc. roy. Soc. A, 231, 211. 

JEFFREYS, H., and Brann, M. E. M., 1951, Mon. Not. R. astr. Soc. geophys. 
Suppl., 6, 148. 

Lust, R., and Scattirer, A., 1954, Zs. f. Ap., 34, 263. 

Nakacawa, Y., 1957, Proc. roy. Soc. A, 240, 108. 

Wastutynski, J., 1946, Astrophys. Norveg, 4. 


[ 859 ] 


The Etching of Diamonds by Low Pressure Oxygen} 


By M. Omar and M. Kenawi 
Physics Department, A’in Shams University, Cairo 


[Received May 9, 1956] 


ABSTRACT 

Diamonds have been etched by heating on a molybdenum filament in 
a low pressure (about 2-5. Hg) of oxygen. The etch figures obtained by 
this method have the same orientation as those obtained by etching in 
oxidizing melts, but differ in outline by being sharply triangular. They 
are in fact similar in all respects, except for the inverted orientation, to 
the natural features called ‘ trigons’. It is therefore concluded, contrary 
to recently accepted opinion, that the latter are also probably etch figures, 
produced by an unknown natural etching process. 


Note: Unless the edge of the crystal is in the field of view, all upright 
triangles are etch pits. 


$1. InrRopUCTION 
Ercu figures on diamond octahedral faces are generally produced by 
heating the crystal in fused potassium nitrate or sodium carbonate at 
temperatures ranging from 525°c to 900°C (see for example Fersmann 
and Goldschmidt 1911, Williams 1932, Omar eé al. 1954). The fuses 
give a limited supply of oxygen which attacks the carbon atoms which 
are diamond. 

The figures obtained are triangular depressions mainly of a pyramidical 
nature with rounded corners, and are oppositely oriented to the natural 
figures called ‘trigons’. Such trigons exist on nearly all octahedral 
faces of diamond, and are a characteristic feature of them. ‘The trigons 
were previously thought to be etch figures but it is mainly due to the work 
of Tolansky (1947 onwards) that trigons are now generally accepted as 
a growth phenomena. The present work is however a retreat back 
to the old view that trigons are most probably etch figures. 

Tolansky was not the only worker who attributed trigons to growth, 
but was preceded by illustrious figures such as Crookes (1909), Fersmann 
and Goldschmidt (1911), Van Der Veen (1911) and Sutton (1928). To 
confirm their view, Tolansky and Wilcock (1947) used multiple beam 
interferometry under what they called ‘high dispersion’ technique. 
They discovered certain marks which appear as extensions to the sides of 
trigons. These extensions, called by Tolansky ‘ basal-extensions ’, 
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were found to be composed of groups of steps leading partly to the 
trigons. The steps were considered as growth layer edges, and as such 
were taken as a definite sign of growth. ‘Tolansky went so far as to incor- 
porate these extensions in a theory, which according to him explained 
the reason of formation of trigons. Williams (1932) wrote at length 
supporting the growth view from a different angle. He presented series 
of photomicrographs which showed that certain hillocks abide by most 
of the trigons. Generally 2 or 3 hillocks confined a trigon. The hillocks 
were composed of layers and were described by him as ‘ growth hillocks ’*. 
He deduced that trigons must be attributed to the same cause, i.e. growth. 


§ 2, OBSERVATIONS 
Experimental procedure 


Both characteristic features of assumed growth, i.e. “ basal extensions ’ 
and ‘ growth hillocks ’ have been produced artificially in a series of experi- 
ments which we have carried out on the etching of diamond using a new 
technique. For this purpose the crystal was placed on a molybdenum 
filament under vacuum of the silvering plant and heated to a temperature 
of about 1000°c. The favourable pressure used was about 2-5 microns 
of mercury. After a long time of etching (about seven hours) etch figures 
made their appearance. The figures obtained got gradually bigger and 
bigger, and smaller figures appeared, until finally a pattern was obtained 
which could not be differentiated from the usual pattern of trigons. This 
observation is new and has not been described before. 


The New Etch Figures 

Figure lj is such a pattern obtained after 21 hours of etching. 
Figure 2 shows the orientation of these pits with respect to the face 
edge of the crystal. As expected the corners of the pits point towards 
the corners of the triangular crystal face. Figure 3 shows a group of 
natural trigons, and fig. 4 their orientation with respect to the crystal 
edge. As can be seen from figs. 2 and 4 the orientation of both sets of 
triangles is different while the grouping of both (cf. figs. 1 and 3) is 
identical. Both types of triangles are straight edged with no rounding at 
the corners of the etch pits as is indeed clear in figs. 1 and 2. 

The above magnifications are not enough to show the two characteristic 
features referred to above, i.e. basal extensions and ‘ growth’ hillocks. 
A, B, C and D, fig. 5 are the basal extensions discovered by Tolansky in 
the case of trigons, while E, F, fig. 6 is a similar extension which occurred 
to an ‘ etch pit ’ while the crystal was under etching in an atmosphere of 
oxygen under reduced pressure. The basal extension of fig. 6 appears 
as a sharp-edged shelf exactly as has been described by Tolansky in the 
case of trigons. Sometimes it extended as a continuation to two sides of 
a single etch pit, as is clearly shown, also in fig. 6. This latter observation 
is a characteristic feature of most of the trigons. 


rk ee 
{ All figures are plates. 


Etching of Diamonds by Low Pressure Oxygen 861 


It has been mentioned that diamond octahedral faces are characterized 
not only by the existence of trigons, but also by the frequent occurrence 
of ‘ growth’ hillocks by the sides of the trigons. These are shown as 
A, B, and C, fig. 7 which is a feature of natural diamond. For the sake 
of comparison we present fig. 8 which has been chosen from our series of 
photomicrographs of etching. X, Y and Z, fig. 8, are hillocks surrounding 
the etch pits in a manner similar to those in fig. 7. While A, B and CG 
appear naturally on octahedral faces of diamond. X, Y and Z are the hillocks 
created in the etching process ; in other words they are ‘ etch hillocks ’. 
The unexpected exact similarity of both kinds of hillocks is surprising. 
Both are composed of layers which have curved edges at the top, and 
except for the fact that they are oppositely oriented with respect to each 
other, one can say that they are similar in all essential respects. 

In view of the above observations, it can safely be stated that, apart 
from the orientation referred to, no morphological difference exists whatso- 
ever between an assumed growth pattern on a diamond crystal and an 
etch pattern created in the manner described above. It is natural to 
attribute both phenomena to the same origin, i.e. to an etching mechanism. 
The side walls of the triangular depressions on (111) are in both cases of 
the form (h k k), (k h k), k k h), with k<k in the natural trigons and 
h>k in those produced by etching in low-pressure oxygen. Regarding 
the reason for this difference no more can be said than that it is well 
known that different etching agents can have different etching habits. 
We know nothing of the conditions under which the diamond has been 
formed, but on account of the well known affinity between oxygen and hot 
carbon atoms, it is relieving to deduce that the etchant in the case of 
diamond cannot have been oxygen. 

Two new characterizing features have been quite recently laid stress 
upon by Tolansky (1955). The first is supposed to be a property of the 
etch figure. It states that etch pits are generally rounded ; the second 
belongs to ‘ growth’ trigons and describes their ability to lie in linear 
formations. 

When we etched our crystals in an atmosphere of oxygen under reduced 
pressure, we obtained straight edged triangles, as is indeed clear in the 
photomicrographs presented here. It is true that when the crystals are 
heated in fusing melts (Williams 1932, Tolansky 1955) or exposed to a 
hot oxygen jet (Shultink e¢ al. 1954) rounded etch figures are generally 
obtained. The rounding of these figures cannot, however, be taken as a 
crystallographic property of the etch pit. It is most probably due to 
solution currents in the oxygen-giving fuses, as well as to eddy currents 
in the jet. 

The property of ‘ growth ’ trigons to lie in straight line formation at 60° 
is not a new discovery. It has already been described in rare cases by 
Sutton (1928). This rare grouping has, however, occurred in a crystal 
that we have etched. The linear formations in fig. 9 are made of 
small etch pits, not ‘ growth’ trigons, and the pits are seen in lines making 
60° with each other, exactly as in the case of assumed growth. 
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The six-sided pit in fig. 9 enlarged in fig. 10 has originally been a trigon, 
inverted of course to the now existing etch pits. The etching extended its 
area, causing a slight truncation at the corners. Further etching caused 
the truncated corners to extend, forming steps simulating growth, while 
the originally straight edges of the trigon are now the curved sides of the 
figure. Ultimately the ‘ growth’ trigon will be converted into a larger 
etch pit, oppositely oriented to the original figure. The importance of 
fig. 10 is that steps on the face of a crystal are not necessarily a sign of 
natural growth. They arise as well during a process of etching. 


The Role of Basal Extensions 

Figures 11 and 12 show a group of etch pits at two successive 
stages of the etching process. They are produced at a somewhat large 
magnification in order to illustrate a very interesting property of the basal 
extension. A,B, C and D, fig. 11, are two basal extensions to two neigh- 
bouring etch pits. Further etching increased their length somewhat as 
in fig. 12. Ultimately they will close, forming a larger etch pit. The 
final pit will contain two smaller pits inside it. The mystery of formation 
of pits inside the territory of larger ones is at last revealed. The basal 
extension turns out to be the dynamic agent responsible for the occurrence 
of this interesting observation, far from being a mechanism of growth as 
has been assumed in the case of natural trigons. 

On account of the unexpected exact similarity of etch and ‘ growth” 
patterns, it is sometimes difficult to discriminate one from the other, 
unless of course the edge of the octahedron is in the field of view. It is 
interesting to study series of photomicrographs of the same etch setting to 
see how the etching proceeds. Figures 13, 14, 15 and 16 illustrate 
different stages in the etching process at 1000°c in an atmosphere of 
oxygen at a pressure of 2-5 microns of mercury. The etching periods 
were 10, 15, 21 and 25 hours respectively. 


§ 3. CONCLUSION 

It is thus clear that every property characterizing the so-called trigons 
or trigon formations has its prototype in etch phenomena. These proper- 
ties are (a) extreme linearity of the edges of the triangular figures, (b) 
the step formation of the walls of the pyramidical depressions, (c) the 
occurrence of basal extensions, (d) the existence of hillocks by the sides 
of the triangles, (¢) the ability of the triangles to lie occasionally in linear 
formation at 60°, (f) the existence of smaller figures inside the larger ones. 
It is difficult to attribute exactly similar features and formations, 
especially if they are of an odd nature, to different causes, and we conclude 
by attributing all triangular depressions and ‘ apparent’ elevations in 
the form of hillocks by the sides of the ‘ trigons ’ to an etching mechanism. 
The word apparent is used because the hillocks abiding by the trigons are 
not actual elevations, but it is what is left over after the etching proceeded 
downwards forming the trigons. Since artificial etch pits, using oxygen 


Etching of Diamonds by Low Pressure Oxygen 863 


as an etchant at a moderately high temperature, are inverted with respect 
to the trigons (which in view of the authors are natural etch pits) we 
conclude that the natural etchant could not have been oxygen. Oxygen 
figures—oppositely oriented to the natural figures—are nearly non- 
existent on the faces of natural diamonds. Williams (1932), after sorting 
out hundreds of thousands of diamond octahedral crystals, found some 
such figures on the faces of a single crystal. It is thus clear that diamond 
grew in an environment devoid of oxygen. If the trigons are etch figures, 
an opinion which we strongly support, some other natural reagent is 
responsible for their existence. 
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SUMMARY 


It is shown that when free carriers are generated in a solid containing 
two or more types of defect capable of trapping electrons or holes, the 
occupancy of one type of trap is not, in general, a simple irradiation time 
dependent function of the form 1—exp (—t/r) but is dependent on the 
parameters of the other traps and at saturation can be significantly less 
_than unity. 

Using a model for crystalline quartz it is possible to explain radiation 
bleaching by both x-rays and neutrons, by considering the dependence 
of the occupancy of one trap on other traps. With the same model an 
explanation of the effect of heat treatment on the kinetics of colouration 
observed by Forman (1951) is given. From a quantitative comparison 
between the predicted and observed rate curves of an X-ray irradiated 
crystal, f-values and defect concentrations are calculated. Less detailed 
comments are made on the coloration of fused quartz. 


§ 1. INTRODUCTION 


OBSERVATIONS on changes in the concentration of colour centres have 
become part of the standard method of gaining information about defects 
in certain solids. In some instances the defect forms a colour centre and 
so gives rise to absorption only after it has trapped an electron or hole. 
The concentration of such colour centres, and hence the absorption, will 
depend on both the concentration of defects and the fraction of them 
which have trapped the requisite charge carriers. It is important to 
know, therefore, the factors which determine the fraction of traps which 
are filled ; they are considered in this paper. In particular it is demon- 
strated how the kinetics governing the trapping of a charge carrier at one 
type of defect may be influenced by trapping at another, or several other 
types of defect. 

The analysis has been carried out for quartz though the results obtained 
may be more generally applicable. Quartz has been chosen to simplify 
the analysis since it is known that defects are not created or annealed. nor 
are colour centres thermally bleached during x-ray irradiation at room 
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temperature (Mitchell and Paige 1956). Data relating the change of 
absorption to radiation dose (rate curves) have been published by Forman 
(1951), Ditchburn et al. (1954) and Mitchell and Paige (1956) for crystalline 
quartz and by Dainton and Rowbottom (1954), Levy and Varley (1955), 
and Ditchburn et al. (1954) for fused quartz. In some of the rate curves 
the absorption passes through a maximum before reaching a constant 
value. This unexplained feature has been called radiation bleaching. 
It is shown by the B, band of one sample in fig. 1, which presents some 
previously unpublished rate curves of two specimens of fused quartz 
X-ray irradiated under identical conditions. By considering the depen- 
dence of the occupancy of one trapping state on that of the other trapping 
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states it will be shown that radiation bleaching, together with the effects 
of previous heat treatment on rate curves observed by Forman (1951), 
may be explained. It will also be demonstrated that there is no justi- 
fication for assuming the fraction of occupied trapping states to be unity 
when the absorption saturates. Often this assumption is made implicitly 
in calculating defect concentrations. 


§ 2, A ‘ SINGLE CENTRE ’ MODEL 
To introduce the notation and the method of analysis together with 
some general assumptions, a simple model is considered. In this model 
only electrons can be trapped to form colour centres. Holes are immobile. 
For this and other models used it is assumed that traps are randomly 
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distributed and that phonon induced transitions between the trapping 
state and an energy band are negligible. 

When quartz is irradiated with x-rays we are concerned only with 
electronic processes such as the generation of free carriers, the occupancy 
of existing trapping states and the recombination of electrons and holes. 
If Z is the number of photons incident on unit area in unit time, »., the 
photoelectric absorption coefficient and y the number of electron-hole 
pairs created by each photoelectron then the rate of generation of electron- 
hole pairs, X, is given by Iu,y for specimens of thickness 2 such that 
/t»- <1. For such samples X can be assumed uniform throughout the 
volume of the crystal. Wolff (1954) has shown that an upper limit to 
y is 2H/Eg where £ is the photoelectron energy, and Hg the gap width 
of the solid. For 30kv x-rays and a typical dose rate X is about 
21014 em sec-!?. The free electron concentration is reduced by 
recombination and trapping processes. The trapped electron can be 
released by either a sufficiently energetic electron or by direct interaction 
with the incident radiation. Release by both these mechanisms is 
assumed negligible. The differential equations governing the free electron, 
trapped electron and hole concentrations are 


dn/dt=X—C n(No—no)—knp, . . .. . (1a) 
dng /dt=C n(N o—N¢), it ya eae te ee 
apidt=X—knps —  . Se 
‘The consideration of charge neutrality implies that 
a Nog=P- a a, Fale 


Here n, p, n¢ and NV, are the concentrations of free electrons, holes, trapped 
electrons and electron traps respectively. C, and k are rate constants given 
by (vo ,) and (vc) where v is the electron velocity, o, is the cross section 
for the capture of an electron at a defect and o the recombination cross 
section of an electron and hole. The v, o product is averaged over the 
free electron distribution. The choice of symbols in this and the next two 
sections has been made in anticipation of their application to quartz. 

Tn order to simplify the solution of these equations, certain assumptions 
will be introduced which may be qualitatively justified for quartz. It has 
been observed by Choong (1945) and others that only fluorescence occurs 
when quartz is irradiated with x-rays. Since it is certain that some of the 
fluorescence is due to trapping or recombination of free carriers it can be 
assumed that the free carrier concentration decays in a fraction of a 
minute. Using a typical dose rate coloration takes many hours. Thus 
the free carrier concentration is governed by a very short time constant 
compared with that controlling the formation of centres. Therefore we 
may regard the free carrier concentration to be at an equilibrium value 
‘determined by the occupancy of traps. 

In the absence of traps eqns. (1 a) and (1 c) become identical and may be 
solved directly to give a 


N= p= Xp tanh (t/7), 


Kinetics of Colour Centre Formation in Quartz 867 
where yx, is the saturation concentration of free electrons given by 
[X/h]"? and + a time constant of [k.X]-1/%. For a typical value of 
X=2x10'4em-? sec}, taking v=107 cm sec! and c=10-!8 to 10-22 em?, 
values of y=4xX10"%cm-3, 7r=2X10-2sec and y,=4x10!4 em-3, 
7=2sec are obtained for the respective values of c. We note that this 
time constant is short compared with hours and that values of y, are 
small compared with detectable concentration of centres in thin samples 
(101° cm? or greater in a 1 mm thick crystal). The presence of trapping 
states will lower the equilibrium free electron concentration. Any changes 
therefore, which occur in the population of traps after the x-irradiation 
has ceased are insignificant. 

In the presence of V, electron traps the equilibrium concentration of 


free electrons is 
Ay,2\ 1/2 


_ OW etncelk—C)) 
i as 
Inspection of b shows that when 7, reaches an observable concentration, 
b>yxo, so that y=y_2/b. Substituting this value of y in eqn. (16) and 
integrating, using the condition that n,=0 at t=0, gives 

In (l—y)+y(kK—C)/k= —XCjt/(N¢ - k), 
where y is the fraction of occupied traps. If k=C,, then the solution 


reduces to 
(-om[-£.) : 
y= 1—exp. | — —— .t-)]: See. Mee Re 
y p No 


This is the simplest rate law which can govern the formation of colour 
centres. However, it is seen here to occur as a special solution for a very 
simple model. The time constant of eqn. (2), N,/X, is 1x 104 sec for a 
trap concentration of 2 101% cm~? and the value of X previously chosen. 
This exceeds the time constant governing the free electron concentration 
by several orders of magnitude, an essential result for the self consistency 
of the method of analysis. 


where 


b 


§ 3. Two TRAPPING STATES 


Two trapping states are present in this model, one traps an electron, the 
other a hole. Both charge carriers are mobile. The main purpose of this 
model is to demonstrate how the parameters controlling the trapping at 
one type of defect may affect the kinetics governing the trapping at 
another. This interaction is transmitted by the free electron and hole 
concentrations, both of which may be controlled by the presence of 
trapping states. 

Recombination of electrons and holes can take place in three ways ; by 
(lirect recombination as in the previous model, by a free electron recombin- 
ing with a trapped hole (rate constant 4,) and by a free hole recombining 
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with a trapped electron (rate constant C,). The differential equations 
for the concentration of centres and free carriers are 
dnjdt=X —C n(No—nc)—A gh 4—knp, ( 
dp[dt=X —A,p(P 4—P4)—C oPNc—kmp, ( 
dno|dti=C nN p—Ne)—C apne, SS nek leas enon oaCy 
dp lat=A,p(P4—p4)—A mpg ( 
To conserve change neutrality, assuming the defects are charged, 
Ten == ppg VS ee eee one 
IN eed 2 
P., is the density of hole traps and p, the concentration of holes trapped 
at them. <A, is the rate constant for the trapping of a hole at a vacant 
hole trap. 

Terms representing the excitation of electrons and holes from traps by 
direct absorption of an x-ray photon or by the impact of a high energy 
carrier have been omitted from eqn. (3). So have terms giving the rate 
of formation of centres by similar excitation processes occurring in the 
vicinity of a defect. 

If C,=C,=A,=A,z, then for any value of & the equilibrium concentra- 
tions of free electrons, y, and holes, %, are independent of the fraction of 
filled traps, 

x=P= A/C No. 
The values of y and % are obtained by the same method of approximation 
used in the previous section. It may be shown that this approximation 
is valid provided C,?, O,?, A,?, 4,2 10-*8k, a condition which can be 
expected to be fulfilled in a real crystal. 

The equations for n, and p, reduce to the simple expression, 

No=Ps=tN o{1— exp (—2Xt/N,)}. . 2. 1h 2. 4) 
If however C,, Cy, A, and A, are not all equal the values of y and 7 are 


functions of the concentrations of centres formed. Suppose C,=C,=4,=hk, 
but that 4,>h, then the equilibrium free carrier concentrations are 


X=A/[KNo+(A2—k) pq], 
AH MAN ce 
Substituting these values in eqn. (3d) and integrating, using the condition 
that p,=0 at t=0, gives 
pa={Y(L+h/[4.—k])—¥ exp (—2Xi/N)}2—YINg, . (5) 
Y is kN,2/[A.—A]. 
The saturation density of centres is 
Nol+ V(A,/-)]. 


Similarly, when C,=A,, C,=A, and O,AC,, a more complicated rate 
law than eqn. (4) is obtained and the saturation concentration of centres is 


No=P4a=No/[1+C,/Cy]. 


where 


Kinetics of Colour Centre Formation in Quartz 869 


From these examples we see that an essential condition for a simple rate 
law to be followed is that the free carrier concentration must be inde- 
pendent of the fraction of occupied states. This is achieved only when all 
the rate constants for transitions associated with the traps are equal, a 
condition which is unlikely to occur in a real crystal. When the free 
carrier concentration is sensitive to the density of colour centres the rate 
law is complicated and the saturation density of centres becomes a 
function of the rate constants. Only if the release of trapped electrons 
and holes is unimportant (values of A, and C, very small compared with 
A, and C,) can the saturation concentration of centres equal the density 
of defects. 


§ 4. A MopEL FoR CRYSTALLINE QUARTZ 


Three absorption maxima have been observed in x-ray irradiated 
crystalline quartz by Mitchell and Paige (1956). A more realistic model 
for quartz, therefore, is obtained by introducing another trapping state, 
of density Pp. into the ‘two centre’ model previously considered. If 
these defects trap holes, the rate of colour centre formation is 


dppl[at=D,p(Pp—Pp)—DenPp: fe tet 6G) 


where p, is the density of D centres and D, and D, are rate constants 
corresponding to A, and A, for the A centre. Equations (3c) and (3 d) 
are unaltered, while (3 a), (3b) and (3 e) become 


dn[dt=X —C n(N o—No)—n(AapyztDepp)—knp, -~ . . (66) 

dp[dti=X —A,p(P4—p4)—D w(Pp—Pp)—Copng—knp, . (6c) 
No=P ath, 

Wa Noe=PI-ParPp. 2. «oe « « « ~» (6d) 


Inspection of these equations shows that the same general conclusion for 
the two-centre model is valid, namely that only when the free carrier 
concentration is independent of the density of centres are the rate laws for 
all the centres of the simple form given by eqn. (4). If the rate constants 
are not all equal the kinetics of colour centre formation become more 
complicated. To consider an interesting example, suppose all the rate 
constants have the same value except D, which is sufficiently small 
compared with the other rate constant for it to be taken as zero. Thus 
holes trapped at D defects cannot be released so the saturation concentra- 
tion of D centres will inevitably equal the density of D defects. If the 
free carrier concentration was independent of the colour centre concentra- 
tion, the saturation values of n,,and p, would be $V, and }P , respectively. 
Clearly, when the value of p, exceeds }P, the influence of the D centre 
will become marked, raising the concentration of C centres above NV, and 
depressing the density of A centres below }P,. It is proposed that this 
is the type of process which gives rise to radiation bleaching. Taking the 
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values of the rate constants chosen above, the equilibrium concentrations 
of electrons and holes are given by 
x d x 
=—- an ee 
Ves ENE X~ k(No—Dp) 


The concentration of D centres may readily be shown to be 


x 
po=Pp(1—exp| —x.-4]). a. Geeectigge Beamon 


Substituting in eqn. (3 d) for x, and pp gives 
dps XX f — (14+2a+2) 
di eh oar oa a 
where z=exp (—Xt/N,), and z<=P,/Pp. Putting «=1 and integrating, 
using the condition that p,=0 at t=0, gives 
p4a=P ,(1+32+322—723)/3(1+2)?, . . . . - (9) 


(8) 


Fig. 2 
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The variation of the occupancy of A, C and D defects with x-ray dose computed 
from eqns. (7) and (9). 


The rate curve represented by this equation passes through a maximum 
before attaining an equilibrium value of }P,. Such a dependence of 
absorption on radiation dose has been called radiation bleaching. The 
variation of the fraction of occupied traps with dose is shown in fig. 2 
for No, DP 4 and Pp Where n¢ was obtained from np=p,+pp. A quantita- 
tive comparison between the observed rate curves and eqn. (9) is presented 
in § 5. 
§ 5. OBSERVED RatTE CuRVES 
5.1. Radiation Bleaching by X-rays 


The three absorption maxima which have been observed in X-ray 
irradiated crystalline quartz have been called the 4, C and D bands. The 
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rate curves for the A and C bands (Ditchburn et al. 1954) show that the 4 
band is radiation bleached. On the present theory, for radiation bleaching 
to occur there must be two defects which trap the same charge carrier. 
Further, it is changes in the occupancy of one of these two types of defect 
which are observed as radiation bleaching. Mitchell and Paige (1956) 
have presented evidence to show that a hole is trapped to form the A 
centre and an electron to form the C centre. Therefore the D defect 
traps a hole to form a colour centre. Because of this conclusion, the A, 
© and PD centres considered in the model used in § 4 are identified as the 
centres giving rise to the A, C and D absorption bands. 

To illustrate the measure of agreement, and the information that can be 
obtained, a comparison between the predicted and observed rate curves will 
bemade. The rate curves for the A and C bands of 100 kv x-ray irradiated 
quartz have been reported by Ditchburn et al.t (1954). The number 
of 100 kv photons incident on unit area per réntgen is 4 10° cm-2R-! 
(Spoull 1946) and their photoelectric absorption coefficient in quartz 


Initial slope in Saturation 
cma? F* absorption in em~t 
observed ox 10-? 4-03 
A band — 
predicted | XP4f410-19/(No. 4H 4) | P4?f410-18/(P4+Nc.4E 4) 
observed feeh0-8 11-1 
C band = 
predicted Xfco10-18/ AKG No. fol0-18/(P4+NcAEKc) 


is 0-5 em~!. Since the crystal used was | mm thick, it can be assumed 
that the rate of generation of free carriers was uniform throughout its 
volume and had a maximum value of 1:7x10%cm?R-!. If AF is 
the half-width of the band (1-8 ev for A, 1-4 ev for C) and pu is the absorp- 
tion coefficient at its peak then p.d#=n. fx 101% where fis the oscillator 
strength for the transition. Using this relation various features of the 
rate curves can be compared with eqns. (3c) and (3d) as shown in 
the table. Equating the observed to predicted values tabulated it is 
foundeethat f,=—1-:0%10-7, fg=5:7K10- and P,—2-1x101? om”, 
No=41x 1019 em, Pp=2-0x 101° cm-%. These are dependent on the 
value of X in such a way that if X is a maximum, minimum values of f and 
maximum values of defect concentrations are obtained. Since P,/Pp=1, 
eqns. (8) and (9) may be used in conjunction with the deduced f-values 
and defect concentrations to predict the maximum absorption of the A 


+ See fig. 12 of this paper. The initial increase of dujdt with dose is due 
to inaccuracies in estimating the dose between the first three measurements 
of absorption. Accordingly the origin has been taken at 210° R in the 
calculations. 
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band and the dose at which it occurs, both can be found independently 
of X. The calculated values are 4:6 cm~}! and 2-9 10° R respectively. 
These are in good agreement with the observed values of 4:35 cm™* 
and between 2 and 3x 10® R. 

Mitchell and Paige (1956) have shown that the A defects are sub- 
stitutional aluminium atoms and that for samples having a maximum A 
band absorption of 5 cm~1, less than 50% of the total Al present in the 
crystal occurs substitutionally. Spectrographic analysis of the total Al 
content showed concentrations of 5 10!®cm-% or more Al atoms for 
similar crystals (Ditchburn et al. 1954). An estimated concentration of 
2x 10!%em-3 A defects from the rate curves is consistent with these data. 

A more complicated model is required for fused quartz since five 
absorption bands have been observed (Paige 1955). However, radiation 
bleaching is believed to originate in the same way as in crystalline quartz. 
Figure | shows that radiation bleaching is dependent on the sample ; the 
bleaching of the B, band is observed for one specimen but not for another 
although they were irradiated under identical conditions. This can be 
explained by a difference in the relative concentrations of various defects. 
In fact, if in eqn. (8), « is made large or small compared with unity, 
eqn. (9) does not pass through a maximum. 

If a simple rate law eqn. (4) is obeyed then a plot of In (u,—) against 
dose is linear where p, is the saturation value of the absorption coefficient. 
Figure 3 illustrates that a simple rate law is not followed by the A 
band of the silica sample shown in fig. (16). <A plot of In (1—(/u,)?) 
against dose is more nearly linear. This function was chosen because if 
A,>k in eqn. (5), then In [1—(y/p,)”] is proportional to dose. 

Rate curves similar to fig. 1 have been observed by Dainton and Row- 
bottom (1954) except that their irradiation was not prolonged sufficiently 
for either saturation or radiation bleaching to be observed. It will be 
seen that the plot of In (u,—j) against dose could be interpreted as the sum 
of two absorption coefficients each obeying a simple rate law. Dainton 
and Rowbottom suggested that both the A and B, bands did follow a 
simple rate law (eqn. (4)) but that they overlapped each other. This 
explanation requires that a high proportion of the absorption at the A 
band peak (2-3 ev) is due to the B, band. However, Yokota (1952) 
has prepared fused quartz which shows a B, band but no absorption 
for photon energies less than 3-0 ev. We may conclude, therefore, that the 
A band is following a more complicated rate law than that suggested by 
Dainton and Rowbottom and that it can be closely approximated by 
eqn. (5). 


J, 


5.2. Forman’s Results 


Forman (1951) has published rate curves for the A band of crystalline 
quartz. After irradiating his samples till the transmission remained 
constant to within 1% (° saturation ’) he thermally bleached them. On 
re-irradiation it was found that the rate of formation of 4 centres was 
faster and that the change depended on the time and temperature of the 
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heat treatment. No change of the saturation absorption and no radiation 
bleaching were observed. Similar heat treatment prior to the initial 
irradiation did not affect the rate of colour centre formation. 

Radiation bleaching may not have occurred in Forman’s samples 
because the relative defect concentrations were not suitable. Alterna- 
tively the condition for ‘saturation’ may have been such that the 
irradiation was terminated before radiation bleaching occurred. 


Fig. 3 


DOSE IN ARBITRARY UNITS. 


Plots of In (us—p) and In (1—[p/p.|*) against x-ray dose for the 4 band of fused 
quartz. 


Since heat treatment prior to irradiation was ineffective, the effect 
observed is electronic in origin; it may be explained as follows. At 
400°C, the temperature at which the heat treatments were carried out. 
the release of trapped carriers from A defects is much more rapid than 
from D defects (Paige 1955). Thus the A band may be completely 
bleached while C and D centres still exist in the crystal. The distri- 
bution of electrons and holes in traps is then completely different from any 
distribution which occurs during the initial x-irradiation. This leads to 
changes in the kinetics of coloration of the 4 band on re-irradiation. To 
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illustrate this, suppose that after a heat treatment all A centres are 
bleached although all D centres are unaffected. With the same condition 
used in explaining radiation bleaching, D,=0, the concentration of elec- 
trons will remain constant during re-irradiation at X .(N~—Pp)* em. 
This is a higher concentration than that encountered during the 
initial irradiation, except during the final stages, hence the A centre 
density will increase more rapidly during the early stages of re-irradiation. 
Longer periods of heating will reduce the D centre concentration signi- 
ficantly and the rate of increase in formation of A centres on re-irradiation 
will decrease accordingly. These are the effects observed by Forman. 


Fig. 4 
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(a) 'The observed rate curves of the A and C bands of neutron bombarded 
crystalline quartz. Concentration of centres in equilibrium with the 
total density of defects. -— — — Before equilibrium is attained. 

(5) Rate curves of the A and C bands of neutron bombarded crystalline quartz 

computed from eqns. (10) and (11). 


Optical bleaching, which can be made more selective than thermal 
bleaching, should be more effective in demonstrating the changes in 
kinetics which can be achieved by disturbing the distribution of trapped 
electrons and holes. 

5.3. Radiation Bleaching by Neutrons 

An alternative means of changing the kinetics of coloration and also the 
equilibrium concentration of centres is to alter the relative densities of 
traps present, or to introduce new trapping states. Mitchell and Paige 
(1956) have demonstrated that fast neutron bombardment of quartz in 
a pile does both ; the concentration of C centres is increased and a new 
absorption band, the # band, appears. The rate curve of the C band of 
neutron bombarded crystalline quartz is shown in fig. 4 (a) with that of 
the A band which is continuously radiation bleached. 
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In the pile the flux of ionizing radiation is such that the rate of creation 
of defects appears to be slow compared with the rate at which electrons 
and holes are trapped. The absorption of a specimen which had received 
a neutron dose of 3-0 101!%n®cm-? was unchanged by subsequent 
X-irradiation. Accordingly the absorption after the initial rapid increase 
is assumed to be due to the equilibrium density of centres determined by 
the total concentration of defects. The continuous radiation bleaching 
of the 4 band can be shown to arise from a change in the occupancy of 
one type of trap due to the creation of other types of traps. The more 
obvious explanation for bleaching by neutron bombardment is that the 
defects in the solid are destroyed (Levy and Varley 1955). It has been 
shown by Mitchell and Paige (1956) that this is incorrect and that the 
bleaching is mainly electronic in origin. 

Suppose that neutron bombardment increases the density of CO defects 
and creates anew defect, #, both at a rate Bt where f is a constant. The 
concentration of A and D defects remain constant. The H defects are 
included in the model for crystalline quartz by assuming that they trap 
holes (to maintain neutrality since the concentration of C defects is 
increasing) and have similar rate constants to the D defects. From eqn. 
(8) the equilibrium concentration of A centres is 

OL Pe 
PACD" 4-1 3a— (PtP BO 
The equilibrium concentration of C centres is 
no(eq)=(Not+Bt}?/(PatNeotft)  . . . . (11) 
Figure 4 (5) show plots of p,(eq) and n,(eq) against dose. A comparison 
of fig. 4 (a) with 4 (b) shows that an explanation of the main features of fig. 
4(a) has been given. Further, when N ,+ft>P,4,n((eq)—Not+ft. 
Thus from the rate of change of 1g with neutron dose for large doses, a 
value of f, can be calculated if a reliable estimate of the neutron damage 
is available. Mitchell and Paige (1956) estimated that 25--5 C defects 
will be produced in quartz per unit thermal neutron dose in a pile. Using 
this estimate and the observed rate of change of absorption an f-value of 
6x 10-3 is obtained for the C band. The variation of no(eq) and p,(eq) 
have been shown as variations of absorption in fig. 4 (5). The absorption 
due to the centres was calculated using f, and f, obtained from the 
x-ray data. 


(10) 


§ 6. CONCLUSIONS 

It has been shown that the occupancy of a trap is dependent on other 
types of traps which may be present. As a result a simple rate law for 
coloration cannot in general be expected. In addition, it is conceivable 
for the saturation concentration of colour centres to differ by orders of 
magnitude from the defect concentration. 

A model for crystalline quartz was set up in which the number of 
types of defects was consistent with the radiation induced absorption 
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spectrum. It was possible to choose the rate constants for the defects 
so as to account for radiation bleaching under x-ray irradiation. Using 
the same model, including the selected rate constants, an explanation of 
Forman’s observations on the effect of heat treatment on the rate curves 
of the A band was given. By extending the model to include defects 
produced by neutrons, the continuous radiation bleaching of the A band 
of neutron bombarded quartz could be explained. From a quantitative 
comparison between the predicted and observed X-ray rate curves, 
f-values and defect concentrations have been calculated. These were in 
agreement with those indicated by radiation damage calculations and the 
known. concentration of aluminium in the crystal. Thus, although a 
major simplification was made by selection of the rate constants and the 
exclusion of certain processes the model is not inconsistent with known 
data. It is important to test its validity further since it provides the 
only means at present of estimating f-values and the fraction of occupied 
defects in quartz. The x-ray rate curves of the D band of a specimen 
which radiation bleached and the effect of thermal and optical bleaching 
on the subsequent rate curves of all three bands would be of great interest. 
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ABSTRACT 


An examination is made of the magnetic behaviour of hematite, 
aFe,O;, in the region of the antiferromagnetic transition temperature 
(—15°c). In this examination it is shown that the magnetic properties 
of hematite may be ascribed to two magnetically independent components 
of ferromagnetism, one dependent upon the antiferromagnetic transition 
whilst the other is independent of this transition. The magnetic proper- 
ties of these two components have been determined firstly from a series 
of measurements with an astatic magnetometer, and then, by direct 
measurement at low temperature (—75°c), the magnetic properties of 
one component have been measured directly. These observations are 
critically discussed with reference to the magnetic structure for hematite 
proposed by Néel, and the inadequacy of this magnetic structure in 
accounting for some of the observations is pointed out. 


$1 

Ty 1950 Morin reported a transition in the magnetic susceptibility of 
hematite powders at about —15°c, and showed that this transition was 
suppressed by the introduction of 1° of titanium ions as an impurity. 
This transition has also been investigated by Néel and Pauthenet 
(1952) who observed the change in the remanent magnetization of a 
single crystal of Elba hematite. The effects on this transition of adding 
small quantities of impurity to hematite in the way described by Morin 
have been more fully investigated and reported in a previous paper 
(Haigh 1957), and these results were compared with the magnetic 
properties of some naturally occurring hematites. In this paper it was 
concluded that the magnetic properties of hematite were not only 
extremely structure sensitive but were also very dependent upon the past 
thermal and magnetic history of the specimens. Also briefly mentioned 
was a cyclic hysteresis effect involving temperature and remagnetization 
which has not previously been reported. 

The present paper will describe in more detail this cyclic hysteresis 
effect with particular reference to the ‘ magnetic-structure ’ for hematite 
proposed by Néel (1953), and further experiments which give more 
detailed information about this magnetic structure will also be described. 


+ Communicated by Professor P. M. S. Blackett, F.B.S. 
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§ 2 

In all the experiments the remanent magnetization of the specimens 
was observed, i.e. measurements were always made in zero field. The 
specimens were in the form of circular discs of plaster of Paris, 2-3 cm in 
diameter and 0-5 cm thick, throughout which 0-5 g «Fe,O, was uniformly 
dispersed. The instruments used in measuring the remanent magnetism 
of the specimens were of two types. 

The first used the astatic magnetometer for observing the changes in 
remanent magnetism of a hematite specimen in a thermal chamber 
situated in zero field immediately below the astatic system. The 
temperature of the specimen was lowered by a flow of cooled air passing 
over it, and by rotation of the specimen in the thermal chamber, changes 
in magnetization could be followed throughout the temperature 
range. 

The second method of observation involved the use of a vibration 
instrument conceived by Blackett and originally built by Sutton. This 
instrument has been described in detail elsewhere (Blackett and Sutton 
1956) ; the method of measurement is briefly as follows: the specimen is 
oscillated at 30 c/s along the axis of two pick-up coils connected in 
opposition. The induced e.m.f., which is proportional to the magnetic 
moment of the specimen, is amplified and recorded with a tuned vibration 
galvanometer. The pick-up coils were placed between the poles of an 
electromagnet and the vibrated specimen was enclosed in a thermal 
box which allowed the temperature to be varied between —70°c and 
100°c. This method allows rapid and absolute measurements of the 
magnetic moment of a specimen in any desired field. It is less sensitive 
but much faster than the magnetometer method. 


oul 


Observations with the astatic magnetometer were made on the 
remanent magnetization of a hematite specimen during the cooling- 
heating cycle. The specimen was magnetized diametrally in a field of 
approximately 6000 Oe, cooled to —80°c and then allowed to regain 
room temperaturey. The variation of the initial remanence, Jp jni¢ With 
temperature is shown in fig. 1. 

As expected, the transition just below 0°c was found on cooling—the 
transition effect being completed at about —60°c. However, on re- 
heating back to room temperature, the remanence did not increase back 
to its original value, but only to about two-thirds of Jrinite On 
remagnetizing the specimen at room temperature in the same field as 
previously, the remanence increased to its former value J R,intt- Lhe 
cooling—heating cycle was repeated and the curves of fig. 1 were again 
traversed during the cycle. Here, then, we have a repeatable ‘ hysteresis- 
effect’ involving change of temperature and remagnetization—the 
portion G to A being completed by remagnetization. 


} Cooling rate ~1-C-1-5 deg./min. 


as 
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§ 3.2 

These experimental observations can be partially explained using the 

“magnetic-model ’ of hematite proposed by Néel, originally for the case 
of a single crystal, but equally applicable to a powder specimen. 

Néel described hematite as being primarily an antiferromagnetic 
body in which the direction of antiferromagnetism is parallel to the 
ternary axis below 260°K, and perpendicular to this axis above—this 
has been confirmed by Shull ef al. (1951) using neutron diffraction. 
However, there is also present at all temperatures an isotropic parasitic 
ferromagnetism whose hysteresis properties are independent of direction 
in the crystal. Superimposed on this, and observable above 260°K 
another parasitic ferromagnetism which is essentially anisotropic and 
dependent not only upon direction in the crystal, but also upon the 
temperature, being closely linked with the direction of the fundamental 
antiferromagnetism. Both these components of ferromagnetism dis- 
appear at the Curie point, which is approximately 685°c. 


Fig. 1 
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To interpret the observed changes below room temperature, Néel 
proposed that the spontaneous magnetization of the parasitic ferro- 
magnetism was in a direction parallel to the fundamental antiferro- 
magnetism. Above 260°K this parasitic ferromagnetism lies in the 
basal plane, but below 260°K, as the direction of antiferromagnetism 
coincides with the ternary axis, so the direction of this anisotropic 
parasitic ferromagnetism is also oriented along the ternary axis. However 


3/0: 2 
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as there is an equal probability of this orientation being in one direction 
as in the other along this axis, the anisotropic parasitic ferromagnetism 
should then disappear. ; 

The cooling part ABCD of fig. 1 is therefore interpreted by Néel in 
the following way: in the region C—D, the observed ferromagnetic 
remanence is due solely to the isotropic parasitic ferromagnetism ; in 
the region A-B, the magnetization is due not only to this isotropic 
ferromagnetism but also to the anisotropic ferromagnetism. In C—D 
the latter averages to zero, in A-B it takes its full spontaneous value. 
The region B-C of the cooling curve may be termed the ‘ transition- 
region ’, in which the parasitic ferromagnetism turns from a direction 
in the basal plane to be along the ternary axis perpendicular to this 
plane on cooling. 


cag 

In the reheating portion of this curve, the main transition of the 
anisotropic component takes place between E and F. As the specimen 
was situated in nearly zero field during this reheating, and assuming for 
the moment that there is no magnetic interaction or coupling between 
the two components of ferromagnetism, then there is no ‘ directive-force ’ 
influencing the anisotropic component to return in a preferred direction 
in the basal plane. Consequently this component should return into 
the plane in a completely random manner and there should be no increase 
in magnetization observed as the specimen heats up through the transition 
temperature. The curve should then follow the direction DEG’ shown 
dotted in fig. 1. The portion DEFG of the curve shows that this is 
not so, and there is a ‘ directional’ return of the anisotropic component 
to the basal plane—this direction being that of the initial direction of 
magnetization of the specimen. 

Considering now the possibility of external directive influences being 
present. By means of the vibration instrument described in § 2, it is 
possible to observe these remanent magnetization variations with 
temperature in the presence of an external field—this introducing a 
directive influence during the reheating. If the cooling and heating is 
carried out in a field ~5000 Oe, then the curve corresponding to ABCD 
(fig. 1) is followed during the reheating as well as during the cooling, and 
Néel’s explanation of the phenomena is valid in both parts of the temper- 
ature cycle. If, however, the observations are made in a weak field, say 
~20 Oe, then the behaviour similar to that shown in fig. 1 is observed. 
As the field at the specimen situated below the astatic magnetometer is at 
least two orders of magnitude less than this field of the vibration instru- 
ment, there will be no directive influence introduced by any incomplete 
cancellation of the earth’s field in the former. 

With regard to the possibility of there being an interaction between 
the two ferromagnetic components, this was also investigated with the 
vibration instrument in the following way. A virgin specimen was 
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cooled in approximately zero field to —75°c when it was magnetized in 
about 17000 Oe. The field was then switched off and the specimen 
allowed to heat up in the small residual field of the magnet (~15 Oe). 
No increase in magnetization was observed as the specimen heated up 
through the transition temperature. There would not, therefore, appear 
to be any directive influence introduced by the isotropic component of 
ferromagnetism. 
§ 3.4 

From the nature of these experiments carried out to determine the 
origin of the directive influence present during the reheating in zero-field, 
it is immediately obvious that it must be a form of ‘ memory ’ process— 
the effects only being observed when the specimen has been magnetized. 
before cooling takes place. 

Here, then, is a serious difficulty in the application of Néel’s inter- 
pretation of the magnetic behaviour of hematite in the case when the 
observations of the remanent magnetization are made in zero field. 

The amount of directional return , i.e. G’G in fig. 1 has been found to 
vary considerably with the particular sample of hematite, being least 
for highly pure, homogeneous and well crystalline samples of hematite, 
although so far it has not been possible to correlate this with any other 
property of the sample. ‘Some of the variations found in the magnetic 
properties of different samples of hematite have been reported in an 
earlier paper (Haigh 1957). 

Having observed the effect of cooling on a magnetized specimen of 
hematite, and the behaviour of the two parasitic ferromagnetisms in the 
cooling process, it is now possible to learn more about these two com- 
ponents by magnetically treating the specimen before subjecting it to the 
cooling cycle. 

§ 4.1 


Blackett (1956) has already discussed the effect of placing a * back- 
field’ on a previously magnetically saturated specimen so as to leave 
the specimen weakly magnetized in a direction opposite to that of satura- 
tion, and has shown that the remanent intensity so produced reverses 
when the specimen is heated. By treating a specimen of hematite in a 
similar manner, and cooling it to —75°c, we can derive information 
about the two components of ferromagnetism, for at low temperatures 
we have only the isotropic component, whilst at room temperature, 
when the specimen was initially magnetized, we have both the isotropic 
(J;) and anisotropic (J,) components. It has been shown that J; is 
relatively unaltered during the temperature cycle, and thus the effective 
J,-value at room temperature can be determined. We are therefore 
able to obtain curves of J; and J, against Hp, the demagnetizing field 
field applied to the specimen after saturation. As we are measuring the 
remanent magnetization of the specimen, the coercivity determined 
from these curves is the remanence coercivity, and will be denoted by 
(H¢),. to distinguish from the coercivity H., measured in the field. 
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§ 4.2 


The procedure in each experiment was to saturate the hematite speci- 
men at room temperature and then place it in given reverse fields. The 
percentage of saturation remanence (J,,;,) resulting from this treatment 
and its sign, either positive or negative according to whether the re- 
manence was in the same or opposite direction to saturation, was measured. 
The specimen ‘was then cooled to —80°c, and the variation of the 
intensity observed—a family of curves of this variation of Jinit, after 
varying treatments, is shown in fig. 2. In the table there is given for 
each of the eight curves the demagnetizing field in column 2, and in 
column 3 and 4 the initial remanence, / init, as a percentage of the satura- 
tion remanence (Js), and in e.m.u./g, respectively—these values being 
determined at 20°c. The J;-values in e.m.u./g determined at —75°c 
when the transition is assumed to be complete, are given in column 5. 
The values of J,, calculated from the vector relation Jinjsx=J,-+-J;. are 
given in e.m.u./g in column 6. 


3 Curve Nos. Refer to 
Col I of the Table 


TEMPERATURE °C 


= 50) 


§ 4.3 


The curve plotted from columns 2 and 4 of the table, shown in fig. 3 
Vv, 

is the usual remanent demagnetization curve for ae specimen when 
determined at room temperature. It gives a saturation remanence of 
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0-15 e.m.u./g and a remanence coercivity ~3200 Oe for the hematite 
specimen. (In this case the magnetizing field of 6000 Oe would not 
completely saturate the specimen. True saturation in hematite seems 
to be possible only with fields greater than 20 000 Oe). 

Fig. 3 
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4 | 1425 ~29 +4-4 —3-8 8:2 

5 .| 1913 ~22 +3:°3 —3-9 7:2 

6 | 2400 ~12 +1-9 —3:8 57 

7 | 2963 ~4-0 +0:6 —3-9 4:5 

8 | 3300 ~—1-6 —0-2 —3-9 3°7 


(The sign convention used above is as given in §4.2) Initial saturating 
field 6000 Oc. 


The curves plotted from columns 2, 5 and 6 of the table are the remanence 
demagentization curves for the two components of ferromagnetism in 
hematite. 
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The demagnetization curve for the isotropic component is a fairly 
smooth curve with a saturation value of 5-1x10-?e.m.u./g and a re- 
manence coercivity ~430 Oe. The shape of the curve is that usually 
observed with a material of low coercive force, i.e. magnetically * soft “— 
this component saturating in fields of about 1200 Oe. 

In contrast the demagnetization curve of the anisotropic component 
differs considerably from that for the isotropic component. With 
increasing demagnetizing field, the remanence does not fall-away rapidly 
and linearly from the saturation value of 9-9 x 10-?e.m.u./g but appears 
to be little affected by fields less than 600 Oe. The magnetization effects 
in this low field region will be considered in more detail in § 4.4. For 
demagnetizing fields greater than about 800 Oe the magnetization falls 


Fig. 4 
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off linearly with the field, giving an extrapolated remanence coercivity 
~4800 Oe. Thus this component has a remanence about double that 
of the isotropic component and a remanence coercivity more than ten 
times greater. With this high value of coercivity, the field necessary to 
saturate this component must be very high, probably of the order of 
several tens of thousands of oersteds. 

Comparing the two components of ferromagnetism, the isotropic 
component is magnetically ‘soft’ whilst the anisotropic component is 
relatively ‘hard’. It thus seems probable that it is this predominant 
anisotropic component of ferromagnetism of a magnetically hard nature 
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that prevents us saturating hematite specimens in moderate fields. The 
magnetic properties of hematite therefore will be determined principally 
by the magnetic properties of this anisotropic component. If the 
relative proportions of these two components in hematite samples can 
be varied, more may be learned about the magnetic structure of the 
mineral. 

From neutron diffraction experiments carried out at room temperature 
and liquid nitrogen temperatures, there is evidence of a difference in the 
proportions of J; and J, components presents in samples of hematite 
which show very little decrease in the remanence on cooling and those 
samples showing a considerable decrease. This work however is still in 
progress and will be reported later. 


§ 4.4 

As a result of these measurements the orientation of the anisotropic 
component in the basal plane at normal temperatures was considered 
both in the magnetized and demagnetized states. On account of the 
trigonal symmetry of the rhombohedral crystal of hematite, it was 
suggested that for the anisotropic component a set of three easy directions 
of magnetization, 120°(27/3) apart, may exist within the basal plane. 
The magnetic properties to be expected from such a system have been 
formulated mathematically by Dr. E. P. Wohlfarth (1955) of the Mathe- 
matics Department of Imperial College, to whom I extend grateful 
thanks for helpful discussions on this subject. Briefly the results to be 
expected are as follows: on magnetizing a completely virgin specimen 
to saturation, no remanence should be observed for fields less than a 
certain critical value. This is the field necessary to turn the intrinsic 
magnetization of a single domain through some minimum angle into an 
easy direction of magnetization nearer to the direction of the applied 
field—this minimum angle being a function of the spatial distribution 
of the easy direction of magnetization in the basal plane. Similarly, 
on demagnetizing a saturated specimen with this magnetic arrangement, 
a critical field must be applied before the state of saturation is altered— 
the critical fields in these two cases obviously being equal. 

Results taken from the demagnetization curve of the anisotropic 
component seem to indicate that this critical field is of the order of 800 to 
1000 Oe (AC in fig. 4). Using this value we find that the predicted 
remanence coercivity, from the model on which Dr. Wohlfarth’s calcula- 
tions are made, is at the most 1500 Oe. From curve 2 of fig. 4, the value 
deduced is more than three times this. 

It is relevant to remark here that from the curves of fig. 4 the isotropic 
component is not strictly magnetically ‘soft’, but shows considerable 
hysteresis properties, having a remanence coercivity ~400 to 500 Oe. 
It is, however, soft when compared with the anisotropic component. 

Whilst this model for the anisotropic component may be considered to 
give qualitative results about the shape of the demagnetization curve 
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and leads to a relatively high value for the coercivity, quantitative 
agreement has not been obtained. This model for the anisotropic 
component may require modification but until the variability in the 
properties of hematite samples is more fully understood, little quantitative 
agreement with theory can be expected. 


§ 5.1 

The next sections describe experiments in which the magnetic properties 
of the isotropic component are investigated directly at low temperatures, 
the observations being made with the vibration instrument described 
in §2. With this instrument it is possible to measure directly the 
demagnetization curve for the isotropic component of ferromagnetism in 
hematite and compare it with that deduced in § 4.2 from the magneto- 
meter observations. 

If we cool a virgin specimen to —75°c, magnetically saturate it and 
then demagnetize it at this same temperature, we will be observing only 
the isotropic component. In order to compare this curve with the 
J,-curve of fig. 4, the magnetization was measured in zero field, after 
the demagnetizing field had been applied—this again gives the remanence 
coercivity (H¢)p. 

Fig. 5 
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Figure 5 shows the directly measured demagnetization curve, which 
agrees well with the curve deduced from the magnetometer measurements, 
the coercivities in the two cases being ~530 Oe and ~430 Oe respectively. 
The differences in saturation moments can be attributed mainly to the 
variable properties of the hematite samples, the proportions of J, and 
J, components being different in the two samples—this possibility having 
been suggested already in § 4.3 as a result of neutron diffraction experi- 
ments. This similarity between the two J;-curves can be seen even more 


clearly in fig. 6 where the normalized demagnetization curves have been 
drawn. 
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§ 5.2 

Following the procedure described in § 5.1 for the demagnetization 
curve at —75°c, demagnetization curves were obtained for virgin speci- 
mens of hematite at temperatures of —40, —20, —10, —5, +20 and 
+40°c. The change in shape of these curves as the transition is 
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approached, and passed through, can be seen in fig. as where the may 
points of interest are firstly, in the movement of the knee of the eaty es 
occurring at about 1000 Oe, across the H p axis for the higher temper- 
ature values, and secondly, in the increase in steepness of the curves past 


888 G. Haigh on Observations on the 


the knee. The implications are that a hard second magnetic component 
is increasing in magnitude as the temperature increases from —20°c, 
thus requiring higher reversed fields to produce zero remanence of the 
specimen. The relative hardness of this new component, compared with 
that of the component existing below —20°c, can be inferred from the 
non-saturation appearance of the curves for fields less than 6000 Oe. 
Using these demagnetization curves, figs. 8 (a) and 8 (b) were drawn. 
They show the variation of the remanence coercivity (He), and satura- 
tion remanent magnetization as the temperature increases from —80°c 
through the transition region to +40°c. In fig. 8 (a) the broken curve 
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is drawn so that the maximum of this curve occurs at the same temper 
ature as the maximum of curve 8 (6), although this latter curve is Tether 
uncertain in the region of 0°c. Whilst with different hematite samples 
there may be a slight shift of these curves along the negative temperature 
axis, the general form of the curve is unlikely to vary much. 
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§ 5.3 

The main results of this section of the paper are that the isothermal 
remanent magnetization curve at —75°c obtained by direct experiment 
is in excellent agreement with that deduced from the magnetometer 
observations of § 4.2, and that the method of analysis employed in § 4.2 
of considering separately the two components of ferromagnetism seem to 
be justified in view of this agreement. As was expected from the 
magetometer results, there is a considerable change in coercivity through 
the transition temperature. The conclusion from considerations of the 
relative proportions of isotropic and anisotropic components present 
(from fig. 8 (b)) and of the change in coercivity through the transition 
(fig. 8 (b)) is that the coercivity of the hard (anisotropic) component 
must be of the order of several thousands of oersted, agreeing with the 
value deduced from fig. 4. 

§ 6 

From the observations described in this paper it is clear that there 
are several problems still to be solved before the magnetic properties of 
hematite are completely understood. 

The theory put forward by Néel to account for the magnetic effects 
observed on cooling a hematite crystal has been shown to be applicable 
in the present investigation when hematite powder specimens were 
cooled to low temperatures. However, the theory does not explain the 
increase in magnetization as the specimen regains room temperature. 
The possibility of there being coupling between the two components of 
ferromagnetism has been investigated, but there does not appear to be 
any coupling of a magnetic nature. 

These experiments, whilst contributing to the knowledge of the two 
parasitic components of ferromagnetism in hematite, do not assist 
greatly in establishing the more important question of the origin of the 
parasitic ferromagnetisms. This latter problem must be closely linked 
with that of accounting for the variability of the magnetic properties of 
different hematite samples. 

As in the previous paper, it must be concluded that the results 
embodied in the present paper are more of qualitative than quantitative 
value on account of this variability of the magnetic properties of hematite 
from sample to sample. 


SUMMARY 


In this paper the effects of cooling and reheating on a magnetized 
specimen of powdered hematite have been described and the observa- 
tions have been discussed with particular reference to the two-component 
magnetic structure of hematite proposed by Néel. 

It has also been shown that an analysis of the two components of 
ferromagnetism in hematite can be made in an original way by observa- 
tions at temperatures above and below the antiferromagnetic transition 
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temperature of —15°c. Whilst one component (J,) disappears on 
cooling through this temperature, the other component (Jy) is negligibly 
affected by this cooling so that the effects of any magnetic treatment on 
this latter component at room temperature, can be investigated by 
cooling below this temperature. Having determined the state of the 
J, component by cooling, then by reference to the initial room temperature 
magnetic condition, the magnetic state of the J, component was deduced. 
In this way the remanence demagnetization curves of a powder specimen 
of « ferric oxide have been determined for each of the components of 
ferromagnetism. 

By direct experiment at low temperatures (—75°c) the demagnetization 
curve for the J; component was determined, and good agreement found 
when compared with the curve obtained from the magnetometer measure- 
ments. 

The magnetic properties of hematite were also examined near the 
transition temperature, where the effects of the anisotropic component 
are most marked. 
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Observations on Extensive Air Showers 
It. The Distribution of Charged Particles} 
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SUMMARY 


The distribution of charged particles in air showers containing 3x 10°, 
3x 107 and 108 particles at sea level has been measured and compared with 
theoretical calculations by Nishimura and Kamata. The best value of the 
age parameter S appears to be 1-4. The significance of this result is 
briefly discussed. 


§ 1. INTRODUCTION 


Tue lateral distribution of charged particles in extensive air showers has 
been the subject of many theoretical and experimental investigations. 
The theoretical calculations are concerned with the Coulomb scattering of 
the electrons in a photon-—electron cascade, and the latest calculations of 
Nishimura and Kamata (quoted in Greisen 1956) seem to give good agree- 
ment with experiment. The distribution is obtained in terms of an age 
parameter S, S=1 representing cascades at the maximum of their develop- 
ment, and S=2 old cascades. The electrons observed in air showers must 
arise from cascades initiated by 7° mesons at varying depths, and the mean 
value of S is usually found to lie between 1-2 and 1-4. 

It is important to determine the structure over a wide range of shower 
sizes, since we might expect that changes in the shower development with 
energy will be reflected in the value of S which gives the best agreement. 
As the energy increases, it is expected that the shower will reach its maxi- 
mum lower in the atmosphere, and that the dominant photon-electron 
cascades may have lower values of NS. 

Over the range of distances from the core which are of interest, the 
density varies by four or five powers of ten. Most experiments cannot 
cover this large range of density, and it is convenient to divide the 
experiments into three groups, according to whether the distances investi- 
gated are between 0—2 metres, 2-100 metres, or beyond 100 metres, calling 
the first group experiments on the core of the shower, the second short 
distance experiments, and the third large distance experiments. We are 
here concerned only with large distance experiments. 

Barrett et al. (1952) have investigated the structure at large distances for 
showers containing between 104 and 10° particles at sea level, corresponding 
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to primary energies of about 1014-101% ev. Dobrotin et al. (1953) pays 
reported measurements at 3860 m on showers containing up to 10 
particles, the primary energy being about the same as in the Barrett 
experiment. Barrett ef al. found good agreement with S=1-4-+0-15, 
and Dobrotin found S=1-25-L0-12. In the present experiments we have 
investigated the structure for showers containing 3 x 10°-8 x 10? parley: 
at sea level, corresponding to primary energies of at least 2x 10" ev. 
We discuss the significance of the distributions obtained. 


Fig. 1 
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§ 2, DESCRIPTION OF APPARATUS 


The apparatus has been described in paper I. In all cases showers are 
detected by ‘a’ coincidences, that is by coincidences between the large 
200 cm? counters in the units. The densities are determined from three 
other pieces of information: (i) the triggering of a tray of six 200 em? 
counters placed at A, fig. 1 ; (ii) the triggering of the small ‘ b ’ counters ; 
(iii) the triggering of other * a’ coincidences. (i) is mainly useful at large 
distances from the shower axis, where the density is of the order of 
1 particle/m? ; (ii) at short distances (densities up to 1000 particles/m2) and 
(iii) can be used at intermediate distances. 


§ 3, ANALYSIS OF RESULTS 

3.1. Selection of Showers 
Three shower sizes have been selected for analysis; I, those showers 
discharging just three units ; IT, those showers discharging eight or more 
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units ; ILI, those showers which it is supposed would have triggered 20 or 
more units on an infinite array. In category I, for the points plotted at 
43, 93 and 134 metres, only those showers which discharged three units in a 
small triangle, none of which was on the outside edge of the array, were 
accepted. The densities at these points are obtained from information 
from the *b’ counters which are either in the triangle discharged, or 
immediately surrounding it. For the points marked X in fig. 2, all showers 
discharging three contiguous units were accepted, provided that no unit 
discharged was on the outside edge of the array. For the remaining points 
obtained by the large tray all three-fold showers were accepted, and the 


Fig. 2 


T 


— > 


DENSITY PARTICLES PER SQUARE METRE 
ro 


“ 
T “ 
| 
I 5 : 
x 3 FOLD SHOWERS, “b’ COUNTERS gee ee a ia 
+8 ® TRAY : 5 ss a 
A % @ FOLD SHOWERS, "b’ COUNTERS Eee aot cne : 
fo} Pees ” %a” n ‘. ws, 
° mo» ” RAY SS5 io “ae : 
© % 20 FOLD SHOWERS, “b” COUNTERS io oe ae - 
0-1 ie ae er ee ; > 5, = : a = i = 


a oa T = 4 Ie }. em be ee Oe 


fe he es eee ee el fe Sy fl eee pe 
4550 60 70 8090 100 120 140 160180200 250 300 350400450500 600 700 8009001000 
DISTANCE IN METRES 


axis was assumed to lie at the centroid of the three units. The strict 
criteria for the shorter distances are necessary to ensure that the shower 
axis is reasonably well defined, but they can be relaxed for the greater 
distances. This point will be clearer when the estimates of distance 
are considered later. 

In category II, for the points obtained from ‘a’ and ‘b’ counters, 
showers were accepted if the axis was judged to lie within the dotted area 
on fig. 1. For the points obtained by the tray, all showers were accepted, 
and again, the axis was assumed to lie at the centroid of the discharged 
‘a’ counters. 
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In category III, all showers were accepted provided that enough units 
on the array were discharged to make possible a reasonably precise 
determination of the position of the axis, and an estimate of the number 
of units which would have been triggered on an infinite array. 


3.2. Density Determination 
We assume that the showers have circular symmetry, that is to say the 
density p is a function only of 7, the distance from the shower core. The 
expected fraction, p, of units in a narrow annulus of radius r to be dis- 
charged is then given by 


p,1=1—exp (—pA,) for single counters of area A, (‘b’ counters or 
6 x 200 cm? tray) 


Po={l—exp (—pA.)}® for coincidence counters, area A, (‘ a’ counters). 


From these expressions and the determination of p by counting, the 
density can be determined. In practice, the number of units in a narrow 
annulus is small. Hence we must average over many showers, which will 
not necessarily have the same density at 7, and over relatively wide annuli. 
We first consider the effect of averaging over many showers, and in § 3.4 
we consider the effect of taking wide annuli. 


Table 1 


pA for 
assumed distribution 0:02 O-l 0-2 0-6 1-0 2-0 4-0 


single 
counters 0:02 0-094 0:19 0-51 O-81 1-6 3-0 


pA from |- 
formulae | coincidence 
counters 0:028 0-12 0-21 0-56 0-86 0-16 2-9 


If pA is known to be small over all the showers concerned, we may write 
p=p,/A, or \/(po)/Ag. The first of these is linear, and it is clear that the 
mean value of p which is obtained from the experiment gives p, the mean 
value of p. 

In order to get some idea of the errors involved when pA is not small, 
or for the coincidence counters, we assume some form for the a priori 
distribution of p which approximates the distribution existing in the 
experiment. Then we may calculate p and p for this distribution, and 
compare the values with the values given by the simple formulae. 

We have chosen a distribution of p of the form N(p)dp=kp-*dp, from 
some value py set by the shower selection criteria to infinity. Table 1 
gives the values of pA for the assumed distribution, and calculated from the 
simple formulae. It will be seen that for the single counter measurements, 
no errors greater than 10% are involved if p4<0-25, and for coincidence 
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counters if 0-15<pA<0-6. These limits were observed throughout the 
work except for a few measurements at high density, where corrections of 
the order of 15°% were applied. 

In category I, the distribution of densities is certainly much narrower 
than given by the distribution assumed above, because large showers are 
excluded. Thus no corrections were necessary in this category. 


3.3. Procedure in Reduction of Results 


For the determinations with the tray, concentric circles were drawn 
round the tray as shown in fig. 1. Then during the observation period the 
number of showers in each category falling in the annuli was counted, and 
the number of times the tray was discharged recorded. Using the rela- 
tions above, the density is determined. These results are plotted in fig. 2 
at the mean distance for each annulus. . 

For the determinations at intermediate distances in category I, the 
following procedure was adopted. A number of showers satisfying the 
selection criteria was taken, and concentric circles drawn round the cen- 
troid of the struck counters at radii of 150, 250 and 400 m. The number of 
units lying in these annuli was counted, and the average taken over the 
sample. We have thus obtained the mean number of units, W, lying 
between 150 and 250 metres and 250 and 400 metres, from the axes of 
showers satisfying the selection criteria. The same shower data were used 
to find the rate S at which these showers are recorded by the apparatus. 
A much longer observation period was then used to find the rates C,, CO, at 
which the ‘a’ and ‘ b’ counters were discharged in the same annuli, in 
coincidence with showers of this type. 

Then the fractions p,, p, are given by 


ae oa Ba. 
Pi owt SW 


and the density is determined as above. 


3.4. Estimation of Distance 


The points plotted at 43m, 98m and 138m present the greatest 
difficulty. The point at 43 m is obtained from the number of times ‘b’ 
counters are triggered in the units which define the 3-fold shower. In 
showers in which three units in a small 90m side triangle have been 
triggered, the probability of the ‘b’ counter being simultaneously 
triggered is found to be 0-16, corresponding to a density of 110 particles/m’. 
We have now to determine what distance from the axis this density 
corresponds to. 

- We note that the centre of the triangle is 54 metres from the vertices. 
However, this is the distance corresponding to the measured density only if 
the density falls offas7~1. In fact, rough measurements (by, for instance, 
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plotting these points at the mean distances) show that the slope is about 
1:8. To make a better estimate of distance we proceed as follows. We 
may assume for this calculation that the axes of showers causing the three- 
fold showers of this type fall nearly uniformly over the small triangle. 
(Calculations using an approximate distribution function show that this is 
not quite true, and experimental results on the variation will be given in 
paper V.) Then the probability that a ‘ b’ counter is triggered by showers 
at a distance r from the counter is 


p=[1—exp {—p(")Ay}]- 


The number of shower axes at 7, dr is proportional to r dr. 
Therefore we may write the mean probability 


_ [tex {—plr) Ay] rar 


i= R 
[ rdr 
J0 


where RF is some distance representing the side of the triangle opposite 
the ‘b’ counter. & was put equal to 84 m, but the result is not sensitive 
to the value chosen. 

- By using this expression, and assuming p oc7~18, we find that the point 
110 part/m? should be plotted at about 43 m, instead of 54, and an error of 
about 3m has been allowed for the uncertainties of the calculation. 
Similarly the next point is plotted at 98 m, whereas the distance of the 
counters from the centre of triangle is 108 metres. The correction for the 
138 metre point and more distant points is negligible. 

For categories II and III, the procedure adopted was as follows. For 
each event, the position of the axis was judged visually. Then concentric 
circles were drawn round this point at 50, 100, 200, 300 and 400 metres. 
The number of units lying in each annulus and the number discharged was 
noted. The densities plotted are the result of averaging over 26 showers 
in category IT and 20 in category III. For the densities at greater dis- 
tances, the unshielded tray was used as before. 


3.5. The Density Distribution 
The results obtained are shown in fig. 2. The solid lines represent the 
best parabola through all the points, i.e. the curve is of the form 
p=kr-“thlee, ~The best values for a and 6 are a=0-73, b=0:38. 
This means that at 100 m the distribution behaves like 7-25, and at 400 m 
likes g 24" 


3.6. Zenith Angle Effects 


Up to this point, we have ignored the fact that the showers are distri- 
buted in zenith angle. We show that a small correction can be made for 
this effect, ; 


Observations on Extensive Air Showers : IIL 897 


When a shower axis makes an angle @ with the vertical, particles which 
are at a distance 7,, from the axis measured along the ground are not in 
general r,, from the axis measured normal to the axis, but are at some 
distance between 7,, and 7,, cos 6. To a sufficiently close approximation, 
we can say that the mean r is 4r,, (1+ cos 8). 

If we assume an exponential absorption in the atmosphere, we can write 
the zenith angle distribution as 


J (0) de=J, exp {B(1—sec 6)} sin @ cos 6, dé 


and the average value of r is given by 


7/2 
(1-Lcos 6) (8) dé 
~ 0 

m -a]2 


| 10) a6 


ae £(B)\t 
=in{1+ 5m} ; 


From measurements of the barometer coefficient, we find B~s. Then 
7r=0-957,,. Similarly, the effective area of the cylindrical counters is 
0-95 of the measured area. The main effect of these corrections is to 
subtract 5°% from the apparent total number of particles. 


r=tr 


3.7. The Total Number of Particles 
The dotted curves in fig. 2 are the curves of Nishimura and Kamata for 
S=1-0, 1-4 and 2-0, and 10° particles. The experimental curve appears to 
be in good agreement with S=1-4. If we assume that this distribution 
holds everywhere, the mean numbers of charged particles in the three 
categories are given in table 2. 


Table 2 


Category I II Il 
N 3-210 | 23107 | 9-010’ 


§ 4. Discussion 
The experimental distribution is the distribution of charged particles, 
whereas the calculations refer to electron—photon cascades only. | Measure- 
ments which will be reported later have shown that about 10% of the 
charged particles are j.-mesons (paper IV). However, the distribution of 
u-mesons is so much flatter than that of the electrons, that corrections for 
the -mesons are greater than 10% only at distances beyond 250 m. 


{By(0)= | F exp(—xu) won du. 
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Probably the best determined feature of our distribution for comparison 
with theory is the slope at about 100 metres, i.e. r/ry=1-25 where r, is 
80 m. The theoretical distribution behaves like r~* where « is as 
follows : 

S 1 1-4 2-0 
a 3°2 2°5 1-6 


The observed value is «=2-25-0-4. This may be compared with the 
value of «=2:7 found by Dobrotin et al. (1956) at 3860 metres for 
(i Ny 

With the present apparatus we have no way of determining the 
distribution at distances less than about 50 metres. This is a rather 
important region, because according to the Nishimura and Kamata 
distribution more than one-third of the particles are contained in it. 
However, if the normalization is made at about 70 metres, the determina- 
tion of the mean number of particles assuming the Nishimura and Kamata 
distribution is not sensitive to the choice of S. Moreover, an independent 
method (to be reported later, papers V and VI) of determining the 
number of particles which is substantially unaffected by their distribution 
gives good agreement with the above figures and thus provides confirma- 
tion of the general correctness of the Nishimura and Kamata distribution 
for these showers. 

According to shower theory S determines both the rate of absorption of 
the particles and the structure. It appears that all measurements that 
have been made on the structure and absorption of showers, at all energies 
and altitudes up to 3860 metres, have given the same result, although the 
accuracy of the measurements is admittedly not high. This statement 
applied to showers in the energy region 10% ev has been discussed by 
Greisen (1956). We have now to consider that it is also true for energies 
near 101” ey. 

At the lower energies there are good reasons to believe that the electrons 
detected at sea level come from 7° mesons produced in the nucleonic 
cascade relatively close to the ground (Greisen 1956, Dobrotin et al. 1956). 
It is, of course possible, as some models of the nucleon cascade suggest, 
that at higher energies (>101° ev) the electrons at sea level come from 
7° mesons produced in the first few collisions. It would be possible to 
choose multiplicities in the early collisions so that the shower at sea level 
would have the observed structure and absorption. Against this, 
however, is the fact that at these high energies the absorption would be 
decreasing with increasing energy, in conflict with our measurements of 
barometer coefficients, which are constant, or even increasing with 
increasing energy (paper VII). Moreover, we would expect that the ratio 
of nucleons to electrons would be decreasing with number of electrons if 
the early collisions were making large contributions to the electron 
component at sea level. This is not borne out by the few measurements 
that exist (Nikolsky 1956). 
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ABSTRACT 


The lateral distribution of the penetrating component has been measured 
in showers containing more than 3x 10° charged particles. Between 
300 and 500 metres, the distribution follows a power law with exponent 
1-5, and falls more steeply beyond 500 metres. The median spread of the 
particles is about 500 metres. The interacting component falls off more 
steeply, and was not detectable beyond 150 metres. About 8-5% of all 
the particles are penetrating, and of these, about 10% are interacting. 

The results are compared with theoretical models, and good agreement 
is obtained with the predictions of Rozental, based on a Fermi model of 
the nuclear interaction. The results appear to be in conflict with the 
hypothesis that the primaries of air showers might be high energy dust 
grains. 


§ 1. LyTRopUCTION 


THE penetrating component in extensive air showers is of interest because 
it is more directly descended from the primary interaction than is the 
electron—photon component, which appears at sea level as a superposition 
of many secondary cascades. Theories of the nuclear cascade, which gives 
rise to both components, have been published by Messel (1951, 1954), 
Rozental (1952) and other authors. Among the predictions of the 
theories which may be tested experimentally are the proportion (k) of 
penetrating particles in showers, the numbers of u-mesons and interacting 
particles, and the distribution of these components relative to the shower 
axis. Hxperimental measurements of these quantities have been made by 
many authors (e.g. Janossy and Broadbent 1948, McCusker and Millar 
1951, Cocconi et al. 1949, Kidus et al. 1952). | Most of this work referred to 
showers containing less than 10° particles at sea level. 

The object of the experiment described here was the determination of 
the distribution of .-mesons in showers containing more than 3x 106 
particles at sea level, at distances up to 1 km from the axis; and an 
approximate derivation of the total number of u.-mesons and interacting 


{+ Communicated by the Authors. 
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particles in the showers. The array constructed at Culham by Cranshaw 


and Galbraith (1957, referred to subsequently as I) was used for the 
selection of showers. 


Fig. | 


x4 xt x2 
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X= PENETRATING PARTICLE DETECTOR 


Position of shielded detectors on the array. 


Fig. 2 


. EXPERIMENTAL ARRANGEMENT 


Four shielded ie detectors were used. Their positions on ae 
array are shown in fig. 1, and a diagram of one of the detectors in fig. 2 
It had six 200 cm? G. M. counters, enclosed individually in lead pipes 
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4mm thick. The counters were connected alternately, in two banks of 
three, to quenching and blocking oscillator circuits, and were shielded by 
20 cm of lead above them, 10 cm at the sides, and 5 cm beneath. Methods 
of monitoring the operation of the counters were similar to those described 
in I. 

Either of the single banks could be triggered by a single »-meson of 
energy >450 Mev. Both banks could be struck in coincidence by: 


(a) Two p-mesons. 


(b) An interacting particle producing two or more charged secondaries 
in the lead. 


(c) A single ~-meson producing a knock-on electron or soft shower in 
the lead. 


(d) A single p-meson incident at an angle >40° to the zenith. 


The probability of observing a (c)-type event was reduced to less than 
1% by the lead pipes (Janossy and McCusker 1949) and these events as, 
well as those of type (d), have been neglected in the analysis. The inter- 
action probability in the lead for a nucleon incident vertically was 0-77, 
assuming a geometric mean free path, and the energy required to produce 
two detectable secondaries, of the order of 310° ev (Camerini e¢ al. 
1951). The possibility of electrons or photons triggering the counters 
through a cascade in the lead was neglected, since a particle of about 
10/2 ev is required to produce, on average, one electron at a depth of 
20 cm. 

Signals from the detectors were transmitted to the control unit, and if 
in coincidence with a shower, recorded on the same punched card. The 
shielded units took no part in the selection of showers. 


§ 3. ANALYSIS OF OBSERVATIONS 


Lateral distributions of «-mesons have been derived assuming that : 


(a) Particles are not associated in pairs or groups over an area compar- 
able with that of a shielded detector. 


(0) Showers are single-cored, with horizontal circular symmetry. 


(c) The core of each shower is coincident with the centroid of the units 
struck on the main array, assigning equal weight to each ‘ a’ coincidence, 
and neglecting information from the small ‘ b’ counters. 


(d) Edge effects, arising from the finite size of the array, are negligible. 


The validity of these assumptions will be considered in §5. The 
centroid of each shower was determined with a simple analogue computer, 
which operated directly from the punched cards. 

To determine the distributions, annuli were drawn about each detector, 
and showers grouped according to the annulus in which the centroid fell. 
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Assuming first that u-mesons are the only type of penetrating particle 

present, then for each detector, if : 

N,=number of single bank penetrating events with coincident centroids 
in a given annulus, 

N,=number of double bank penetrating events with coincident centroids 
in a given annulus, 


N,=total number of showers with centroids in the annulus, 


then the density (p) of u-mesons at a distance from the centroid equal to 
the mean radius of the annulus is given approximately by : 


(Ni4N)No={1—exp(—pA)} 2...» () 


where A is the total sensitive area of the detector. Expression (1) 
is not strictly true, since p is distributed within the annulus, and 
{l—exp (—PA)} is not exactly equal to {1—exp(—pA)}, but it has 
been shown previously (Cranshaw et al. 1957, referred to as III) that the 
approximation is satisfactory, provided pA is less than about 0-25. 
For higher densities a correction may be applied. 

Now if only »-mesons were present the number of double bank coin- 
cidences NV, would be given by : 


N JN p—={1—exp (—p4/2)}2 2 . 2... . (2) 


where p is derived from (1). Since interacting particles also occur the 
experimental value of NV, is greater than V,. A better value for the 
jz-meson density is given by : 


(Ny+N,)/Np={1—exp(—pA)} 2. 2. . (8) 


and the mean density of interacting particles, 7, at the same distance from 
the centroid, by : ; 
(N,—N,)/Np={1—exp(—epA)} . . . . . (4) 


where « is the efficiency for detection of an interacting particle as a double 
bank event. Experimentally, both the interacting particle density 7, and 
the correction for the distribution of p, were negligible for distances greater 
than 150 metres from the shower centroid. 

A correction was applied for the chance coincidence rate between the 
shielded detectors and showers observed on the array. This was deter- 
mined from the test cards, produced artificially at half-hour intervals. It 
was 0-25 °%, of the total number of showers observed, for each detector, and 
was subtracted for each annulus, according to the number of shower 
centroids observed within it. The correction to the number of penetrating 
particles was about 5% for most of the annuli, but rose to about 40% at 
distances greater than 800 metres from the detector. 

Possible distortions of the lateral distribution, arising from geomagnetic 
deflections of the »-mesons in flight, were neglected, since the expected 
effects were small ; a study of events from the central detector (3 in fig. 1) 
showed no preferred direction for showers relative to the detector. 
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§ 4. RESULTS 
4.1. .-Mesons 


Approximately 11000 showers have been analysed, with a total of 722 
- penetrating events in coincidence, from all four detectors. The distri- 
bution of these events relative to the detectors have been given in detail 
previously (Porter and Sherwood 1957). Annuli of 84, 105, or 126 metres 
width were used for analysis. Showers were analysed in two groups, 
>3-folds, containing on average 6-7 x 10° charged particles, and >8-folds, 
with 2-7 x 107 particles (see paper III). 

Densities derived from the observations are shown in fig. 3 as a function 
of distance from the axis. The distribution p=Ar‘”+*” has been fitted 
to the points. Expressing p in particles/metre?, and r in metres, the best 
values of K,n and a, are : 

=3-tolds K=612 44 
n=0-825 
a=0-000417 metres” ?. 

Mean Standard Error in (n+ar)=-+0-095. 

The >8-fold results could be fitted to the >3-fold distribution, normalizing 
K, and are statistically compatible with it on a yx? test. The dotted curve 
in fig. 3 is the distribution of all charged particles in >3-folds. 

Total numbers of j.-mesons in the showers can be found by integrating 


r 00 


the experimental distributions: N,= | 27prdr. The distributions are 
“0 


determined only for r=40—1000 metres, and the total numbers of u-mesons 
within these limits, with the contributions added on extrapolating the 


distributions, and the proportion of u-mesons in the showers, are shown in 
the table. 


Number and Proportions of u-Mesons (r in Metres) 


Shower | Experimental ese gad Estimated Fraction of 

size r=40-1000 nn 1000 total number | all particles 

>3-folds | (4:3+-2-1) x 105) 1-5 x 104] 1-3 x 105] (5-7 +2-3) « 105| 0-085 +. 0-035 
(a) 

>8-folds | (1-9+ 1-4) x 10°] 7-6 x 104] 3-3 x 10°} (2-3 £1-3) x 108| 0-085 +.0-049 
(geen O'S) 


The root-mean-square spread of the z-mesons, which is the quantity 
usually predicted by cascade theory, is difficult to measure, since it is 
sensitive to contributions from large distances. In the measured region, 
the r.m.s. spread is about 530 metres, for both >3-folds and >8-folds. 
If the distributions are extrapolated to infinity, the spread becomes about 
1000 metres in both cases. This is probably greater than the true value. 
The median spread may be determined more accurately. It is about 
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500-++ 100 metres for both classes of shower size, assuming the distributions 
to hold beyond r=1000 metres, and about 400 metres in the measured 
region itself. 


4.2. Interacting Particles 


The array is unsuitable for the study of the distribution of interacting 
particles (I.P.) in the showers, since the density falls off rapidly from the 
axis, which is difficult to locate to an accuracy better than about 50 metres. 
The total number of I.P. in the showers may, however, be determined with- 
out deriving the structure function. 

Consider a detector of area A, near the centre of an array of area B. 
If in a time ¢, the number of shower cores detected by the array—C, and 
the number of coincident J.P. detected—=/, then the rate of shower cores 
f.=C/tB per unit area/unit time, and the rate of interacting particles 
R,=I1/tA per unit area/unit time. The rate of cores in an infinitesimal 
annulus of radius 7, dr, about the detector, is dR,=R,x 27r dr per unit 
time. The rate of I.P. at the detector due to showers with cores within 
the annulus, is dR,;,=2z7rR,p dr, per unit area/unit time. The total rate 
of I.P. due to cores anywhere on the array is then 


dele, | Zar para KR, [ 2rrp dr=R,N/unit area/unit time, 
B ~ 0 


where JN is the total number of [.P. in a single shower. The integral over 
the array can be put equal to the integral from 0 to o, since the spread of 
the interacting component is small compared with the array. The total 
number of I.P. in the shower is now N=Rf,/R,=I1B/CA. 

The mean number of interacting particles in >3-fold showers has been 
determined by this method. Only 20 events were available which could 
be ascribed to I.P., from detector 3 in the centre of the array. This 
corresponds to a total number V=(4-1-++-2-1) x 104 particles. No reliable 
estimate could be made for >8-fold events. Further measurements on 
the interacting component are in progress. Interacting particle densities 
could be determined only at two distances from the axis. These are 
shown in fig. 3 with the distribution po7~1->, which is consistent with the 
results. 


§ 5. Discussion 


~ 


5.1. Validity of the Assumptions 

(a) The assumption that penetrating particles are not locally correlated 
is justified by the experimental observation that double bank penetrating 
events occur only in the region where they would be expected from 
considerations of u-meson and interacting particle density, and not in 
parts of the shower distant from the centroid. 

(b) The assumption that the showers are single-cored, appears, from a 
visual examination of the patterns produced, to be valid for the majority 
of events. About 0:5% of the events do not conform with this pattern 
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(Cranshaw 1956) but the effect of including these in the analysis is negli- 
gible. The assumption that showers have horizontal circular symmetry 
is of course invalid, since the zenith angle of incidence may be appreciable. 
The effect of oblique incidence on the observed lateral distribution has 


Density part/m > Bfolds 


Density part m /> 3 fold 


Fig. 3 
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been calculated, using a zenith angle dependence of J=/, cos® 6, for a 
y-meson distribution of the form p= Kr-", with m varying between 0-5 and 
2-0. The zenith angle dependence was determined from observed baro- 
meter coefficients. The contribution from oblique showers leads to an 
underestimate of the true j-meson densities by about 5°, but does not 
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distort the form of the distribution. This correction has not been applied 
to the results given in fig. 3. A more serious source of error may arise 
from a change in the true value of / in oblique showers. The proportion 
of jx-mesons probably increases with zenith angle, since the absorption 
length for .-mesons is greater than that for the electron photon cascade. 
This will give rise to an experimental over-estimate of the absolute number 
of w-mesons in vertical showers, but would not be expected to alter the 
form of the distribution from that which would be observed if only vertical 
showers were present. 


(c) and (d) The assumptions that the centroid of the units struck is a 
good approximation to the shower core, and that edge effects are negligible, 
is justified experimentally by the consistency of the results from different 
detectors. Each of these, if analysed separately, give distributions close 
to the mean, in the >3-fold showers. In the >8-folds, the deviations are 
greater, but since the statistics are poor, the differences may reasonably be 
ascribed to fluctuations. 

The neglect of effects due to knock-on electrons or oblique jz-mesons may 
again be justified experimentally by the small number of double bank 
events observed. Spurious events of this type should occur as a propor- 
tion of the total number of single bank coincidences observed, and since 
very few double bank coincidences occurred for 7 >150 metres, the effect 
should be small also at short distances from the axis. 


5.2. Comparison with Previous Experiments 

The results are in qualitative agreement with previous observations on 
somewhat smaller showers. The increase in the proportion of penetrating 
particles with distance from the axis was reported by Hidus et al., and the 
proportion of interacting particles is consistent with that found by 
Vavilov et al. (1955). 

The form of the »-meson distribution is slightly different from that 
reported by Zatsepin (1956), who found a distribution p=Kr~’, with 
v—0-5 for r=—0—50 metres, and v=2 for r=300—500 metres, in showers of 
105-108 particles. The distribution p=Ar-‘"*“” is tangential at any point 
to a distribution p=Kr—”, where v=d (log p)/d(log r)= —(n+-ar-+-ar In r). 
For >3-fold showers the best value of v for r—300—500 metres is 
v=1-48-+0-06. The slope could not be measured for r—0—50 metres. 


5.3. Comparison with Theory 
The results are in reasonable agreement with the predictions of Rozental 
(1952 loc. cit.) using a model of the nuclear interaction similar to that 
suggested by Fermi (1951, a, b). The ratios of penetrating particles to 
electrons, and of interacting particles to .-mesons, are predicted for showers 
containing 7 x 10° particles at sea level : 
N NV 


ty 
==(()Pille 
N 
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Experimental values for >3-folds are : 


AT N 
sue —0-095-+0-039, —/ —0-07-L0-04. 


AT 
iV 6 =" is 


The comparison of theory with somewhat larger showers is probably 
justified, since only ratios of numbers are considered. The predicted 
r.m.s. spread of the y-mesons is 600 metres, in reasonable agreement with 
experiment. The detailed form of the distribution and its median spread, 
are not given. 

Other models of the nuclear interaction cannot be excluded by the 
experimental results. Unfortunately, detailed predictions of the ;.-meson 
distribution arising from a plural type of cascade (Heitler and Janossy 
1949) or taking into account the tunnel effect in the nucleus (Roesler and 
McCusker 1953) are not available. 

The results seem difficult to reconcile with the suggestion of Alfven and 
Herlofson (1956) that the large showers may be caused by primary dust 
grains, however. In this model the primary contains 103-10° atoms, each 
carrying 1012-10! ev/nucleon. The proportion of p-mesons predicted 
between 70-200 metres is about five times greater than that found 
experimentally, and is only about 50% greater at the core than at 200 
metres. Between 200 and 1000 metres the form of the distribution 
resembles the experimental result, but the proportion of u-mesons is 
again higher. 

The agreement with the theory of Budini and Moliere (1953) is also poor. 


§ 4. CONCLUSION 


The distribution of 4.-mesons in showers containing about 6-7 x 10° 
particles, and produced by primaries of estimated energy 5 x 101° ev has 
been determined with moderate accuracy, between 50 and 1000 metres 
from the shower core. It may be fitted to a law of the form p= Kr-("+™), 
The determination of the total number of u-mesons in the showers, and 
their median spread, is dependent on the extrapolation of the law beyond 
1000 metres, where, if it is assumed to hold, a further contribution of about 
20% is made to both these quantities. The root-mean-square spread is 
sensitive to the contribution from beyond 1000 metres, and cannot be 
reliably determined. The distribution of j.-mesons in the group of larger 
showers, containing about 2-7 x 107 particles, follows a similar law, though 
it may be slightly less steep in the region r=0-400 metres. The propor- 
tion of interacting particles appears similar to that found at lower primary 
energies. 

The qualitative similarity between these results and the work on smaller 
showers seems to indicate that the basic nature of the nuclear interactions 
does not change between 10'® and 5x10! ev. Perhaps the most 
significant feature of the work is the disagreement between experiment and 
the primary dust grain model. The agreement with a theory based on 
Fermi’s model of meson production is reasonably good, but it seems likely 
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that a model in which the interactions are elastic could give equally good 
agreement with experiment. This ambiguity could be at least partially 
resolved if the altitude of ~.-meson production were known. 
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ABSTRACT 


An experiment has recently been reported from which it has been 
inferred that the (p, d) reaction is catalysed by the presence of mu-mesons 
forming mesic atoms and molecules. A brief analysis is made in this note 
of the reaction rates necessary to explain the experimental results. The 
relatively large number of events observed in which the meson takes away 
the energy released in the reaction can be understood as the result of 
conversion of an electric monopole transition. 


$1. [INTRODUCTION 


ALVAREZ et al. (1956) have reported the following experiment. Mu- 
mesons, arising from the decay of pions from the Berkeley cyclotron, were 
slowed down and stopped in a liquid hydrogen bubble chamber. A 
number of events were observed in which the meson came to the end of its 
normal range, and then apparently was re-ejected with an energy of about 
5-4 Mev; on one or two occasions this behaviour was observed twice in 
succession. These events were interpreted as evidence that the nuclear 
reaction 
p+d-— *He- 5-4 Mev 


was occurring with appreciable frequency following the formation of a 
molecule in which the proton and deuteron are bound closely together by 
the meson ; the excess energy was carried off by the meson. 

The experiment was repeated with hydrogen enriched in deuterium ; 
the number of events observed increased with increasing concentration in 
the manner shown in table 1, which contains the data given by Alvarez. 
The first line gives the reciprocal concentration of deuterium, the second 
the number of mesons observed to stop in the chamber, and the third the 
number of events of the above type. The last line is given by a semi- 
empirical formula for the number of events, discussed in § 2. 

When a meson is stopped in condensed material it is captured into the 
Coulomb field of some nucleus, and falls rapidly into the lowest 18 orbit, 


+ Communicated by Dr. B. H. Flowers. Since completion of this paper 
the author has received a preprint of a similar analysis by J. D. Jackson 
(Phys. Rev. (in the press)). 
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as was shown by Fermi and Teller (1947). The mesic atom formed in this 
way with a nucleus of unit charge is a small neutral system that can easily 
diffuse through matter ; the repulsive Coulomb field is largely shielded by 
the meson, so that interactions with other nuclei can occur more readily. 

This shielding effect is not however sufficient by itself to explain the 
observed frequency of reactions. The intrinsic probability of the (p, d) 
reaction should be very nearly equal, at least at thermal energies, to that 
of the mirror (n, d) capture reaction. The latter has a cross section of 
0-46 mb at the neutron velocity of 2-2 x 10° cm/sec which gives 


vo 10-22 em? sec}. 


The particle density of liquid hydrogen is about 4:2 x 102 nuclei per em; 
so that even with complete neglect of the Coulomb barrier the reaction 
rate would only be of the order of a few per second, with a negligible 


Table 1 


I/c 


Mesons stopped 2500 1600 1000 
Events 18 38 30 


Formula (2) 


chance for the reaction during the lifetime of the meson, 2-2 usec. If 
however the proton and deuteron are bound together in a mesic molecule, 
the local density will be increased by something of the order of (207), 
207 being the ratio of mesic to electronic masses; then the reaction can 
occur at a rate comparable with the mesic decay rate. 

These considerations were pointed out by Zeldovich (1954), who also 
made an estimate of the rate at which molecules might be formed. He 
found that the rate was small unless there were a molecular level close 
to the dissociation threshold so that there was a resonance effect, and he 
estimated that the first vibrational level was near this position. A precise 
estimate is clearly not possible without an accurate calculation of the three- 
body molecular system. 

In § 2 we make a phenomenological analysis of the course of the atomic 
and nuclear reactions; we find that the experimental results collected in 
table 1 yield two relations between the various reaction rates. In § 3 we 
consider briefly the atomic reaction rates, and in § 4 the nuclear ones ; 
there the main problem is to understand the frequency with which the 
reaction energy is taken off by the meson, rather than by radiation. 
This can be understood as the result of a monopole transition present in 
addition to the magnetic dipole one that is responsible for most of the 
capture of neutrons by deuterium. 


3224 
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§ 2, PHENOMENOLOGICAL DESCRIPTION 
The lifetime of the mu-meson is 2-2 psec. and the reciprocal, the decay 


rate, will be denoted by 
Apt 0 x 10h see. 

Consider now the history of one meson stopping in the liquid hydrogen, 
with a concentration c of deuterium. After a time negligibly short com- 
pared with the lifetime the meson will be captured into the 1S orbit around 
either a proton or deuteron with probabilities (1—c) and c respectively. 
The binding energies of these mesic atoms are known precisely from the 
usual formula ; there is an appreciable difference between them due to the 
different reduced masses of the meson in the two cases. The binding 
energy of the molecules can only be estimated rather crudely, because the 
adiabatic approximation that is good enough for the electronic case may 
be appreciably in error with mesons. The main part of the molecular 
binding will scale up in proportion to the reduced mass of the meson, but 
the (negative) contribution of the zero-point vibrational energy scales as 
the 2 power of the mass. These estimates of binding energy are collected 
in table 2. 


Table 2 


Mesic atom or molecule-ion 


Binding energy, ev 


The mesic deuterium will end its life either by decay of the meson or by 
formation of a molecular-ion. In addition the mesic hydrogen may ex- 
change its meson with a deuteron, and this is likely to be much faster than 
molecule-formation. The rate for this exchange reaction will be propor- 
tional to the concentration c, and will be denoted by cA,. 

The rate of formation of the different molecules will be different apart 
from the concentration factors because of the effects of the different masses. 
These rates will be denoted as follows : 


for (H), to form (HH),+: (1—c) Ans, 
for (H), to form (HD),+: ca 
for (D), to form (HD),*+: (1—c) A,,3, 
for (D),, to form (DD)*+: 


Hie 


ma: 
Once a mesic molecule-ion has been formed it will persist until either the 
meson decay or a nuclear reaction occur. 

With the reaction rates defined above the number of (HD)+ formed 
from one meson stopping is easily found to be 


Nal C(1—C)A ng [Act Amst 2Ag+(1—C)Ame+CAmal 
ot (ens HelmallAgHoUF Amt Oona]  ° oD 
In this formula we have neglected the possibility that any nuclear reaction 


N 
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may occur before the molecule is formed. Zeldovich (1954) points out 
that the same resonance effect that enhances the rate of molecule formation, 
will also enhance the chance of a reaction ‘in flight’, that we found 
previously to be small ; but the yield of this mechanism is small compared 
with that from the molecule. 

The energy released in the nuclear reaction may be taken off either by a 
gamma ray or by the meson ; the rates for these two processes in a mole- 
cule will be denoted by A, and A, respectively. If the reaction is radiative 
the meson in the molecular orbit will almost certainly be retained in an 
orbit around the final He nucleus, because the recoil velocity of that 
nucleus is small, being only about 1/10 of the orbital velocity of the meson ; 
there it will remain until it decays. On the other hand when the meson is 
ejected with 5-4 mev it will go again through the same cycle of reactions. 
Accordingly the number of events (in which a meson is ejected) to be 
expected for every meson coming into the chamber equals 


Ag tA +A(1—N) 
For small concentrations c, we can write then 
1/Y=A+B/e lew oe Pahhet Poe coli aio g ') | 


where A and B are certain combinations of reaction rates. Provided 
A, is large compared with A, and the 4,,,,, we find 


bi 


A+1=[(Ap+Ams)/Ans] : [(Ag +A +Ao)/Aa] : Heels) 
and B=(A+1)(AgtAme)/Ae- 
The data in table | can be fitted, within the statistical error, by 
A=32, B=0-025 eee eee aC) 


as is indicated by the entries in the last line, obtained by multiplying the 
number of mesons stopping by Y(c). 

Two qualitative conclusions may immediately be drawn from these 
numbers. The first is that \, must be considerably greater than either A, 
or X,,9; this correlates with the onset of a saturation effect in the number of 
events at quite low deuterium concentration. The second is that 

AglAy 2 1A=0-03 
which means that internal conversion is quite large ; equality will hold 
only if both X,,,3, the rate of molecule formation, and 4,, the reaction rate, 
are large compared with the decay rate Ay. 


§ 3. Atomic Reaction Rates 


The second of eqns. (3) implies that A, must be at least of the order of 
109 sec-!; which means that for the reaction 


(H),,--d — (D),,-- p+ 135 ev 


at zero energy the value of vo must be at least 10~' cm see™!. The only 
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simple method of estimating such reaction rates is by the use of Born 
approximation, and in the present case that method leads to the value 
vo=2 10-11 cm? sec-!. It is well known however that the Born 
approximation usually leads to a gross over-estimate of cross sections at 
low energies (see, for example, Dalgarno and Yodav 1953). 

It is also evident from eqn. (3) that however strong internal conversion 
may be, the rate of molecule formation A,,, must be greater than 
\p/(A-+1)=1-5 x 104 sec-1. The energy liberated by the formation is 
most likely to be carried away by an electron, and a reasonable model would 
be provided by the collision of a mesic deuterium atom with an ordinary 
hydrogen atom, 

(D),+(H), > (HD),+-¢+250 ev 

using the rough estimate of binding energy adopted in table 2. An 
estimate of the reaction rate by Born approximation (using a rough 
model of the molecular wave function) gives vo~10-73 cm’ sec}. The 
much greater rate apparently implied by the experiments, must be due 
to transitions to excited molecular states more nearly in resonance with 
the initial state, as was suggested by Zeldovich (1954). We have not 
tried to estimate these more closely. 


§4,. NucLEAR Reaction RATE 


In this section we estimate ,, the rate for the radiative reaction. We 
assume that the intrinsic reaction rate is the same as for the mirror 
reaction, with vo= 10-22 em? sec"! ; the molecular wave function is some 
function (x, 7) where 7 is the separation of the nuclei and w is the distance 
of the meson from their mass-centre. Then 


\=Pue=10-@P |. |... CB) 


where P= f|(x, O) |? d’x in the ‘ collision intensity ’ for the nucleons. 
If further we make the adiabatic assumption 


p(x, r)=(x; 7) x(r) 
where 4(z; 7) is a normalized wave function for the meson with fixed 


nucleons and y(r) is the wave function of the nucleon motion in an 
effective potential V(r), then simply 


P={x(0)[?. se ene a Maca i (Gs) 
Near the position of equilibrium, at the separation 
p=1-064/207=5 x 10-14 em, 


the potential V may be approximated by an oscillator scaled up from the 
electronic case, giving 


V(r)=(207) [—2-8++5-4 (r—p/p)?] ev. 
The normalized ground state wave function is then 
x(r)=1-} (w¥?/4 19? p)1? exp [—hu (r—p/p?] . . (7) 


where the numerical value of w is approximately 3:0. Near r—0, x(r) 
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must behave like a Coulomb wave function; we have joined (7) onto a 
Coulomb function at 7/p=0-36, with a value of 7=1:44; the latter was 
determined by the strict use of the adiabatic approximation. This led to 
a value for P according to eqn. (6), 


P=~10?8 em-3 


and so according to (5) gave A;=108 sec~!, twice the rate of meson decay. 

This estimate has been derived by very approximate analysis, but 
that error is probably less than the error of the adiabatic approximation, 
which should tend to over-estimate the reaction rate. 

It has also been implicitly assumed above that only the S-wave contri- 
butes appreciably to the reaction. There is however a small but signi- 
ficant contribution from the P-wave, because the relatively greater value 
of the electric dipole matrix element partially compensates for the 
reduced penetration of the wave. 


§5. INTERNAL CONVERSION 


The remaining problem is to explain the value of A,, describing the con- 
version process actually observed in the experiment. As we saw in § 2, 
the internal conversion ratio is at least 0-03 from the analysis of the data. 

The chief radiative capture process involves a magnetic dipole transition; 
according to Austern (1952) the amount of mixed electric quadrupole 
(arising from D-wave components of the wave-functions) in the neutron 
capture reaction in deuterium is of the order 10-5, and it cannot be much 
greater for the mirror reaction. However the internal conversion coeffi- 
cient for M1 radiation of 5-4 Mey and for a mesic He atom is only 10-4 
(according to Born approximation), and the quadrupole admixture is too 
small to change this appreciably. 

Thus the important conversion process cannot be that of the dominant 
radiative process. As was mentioned at the end of § 4, there may be a 
certain amount of E1 radiation due to the nuclear P-wave ; since however 
the molecule as a whole is in an S-state the conversion would then be of 
a P-wave meson, and therefore far too small. 

The remaining possibility arising from the nuclear S-wave is conversion 
of the (forbidden) electric monopole radiation. As has recently been 
pointed out by Church and Weneser (1956) this process may compete 
effectively with other radiations in suitable circumstances. In the 
present case, since the final nuclear state is doublet, the nuclear spins 
will be correctly aligned for } of the time. This leads to the following 
expression for the rate of conversion for a meson initially in a 1S He atom 
orbit, according to Born approximation, 


A=(1/3) P (32/9)(e?/fic)?(w/a)|(f |27,7|7) [2 - . - (8) 
where w is the frequency of the ejected meson, a is the mesonic Bohr 


radius, and the monopole matrix element is between the final *He state 
and an initial state appropriate to a p—d collision at zero energy and from 
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which the penetration factor P has been removed ; this initial state can 
therefore be identified with that for n—d scattering at zero energy, nor- 
malized to unit density at large distances. 

We have estimated the monopole matrix element in the following way : 


(a) The final state has been taken in the form const [exp ( —yp)/p] where 
p2 is the sum of squares of distances between particles, and y has been 
selected to give the right asymptotic behaviour for dissociation into proton 
and deuteron, since that is the most important part of the wave function 
in the matrix element. 

(b) The incident wave, at zero energy, is simply the deuteron wave 
function, for which the asymptotic form has been used with the correct 
normalization. The initial state has been put equal to this incident 
wave less a suitable multiple of the final state, so that the two are ortho- 
gonal. 


If the deuteron wave is asymptotic to [exp (—ar)/r], the constant y 
determined in the above way is equal to «/0-6. With these assumptions 
the matrix element can be evaluated analytically, with the result that 
from formulae (5) and (8), 


NA,=0-8 x 10-4 (m, cf)" | (f|X7,2|1) P~O-17. . . . (9) 


This is larger by a factor 6 than the minimum value deduced above ; 
much of this factor is probably due to the over-estimate of the meson 
density at the origin made in the adiabatic approximation. That 
density is about 20 times the density at the same point in the average 
molecular state. 

It is interesting with this value of the monopole matrix element to 
calculate the corresponding cross section for pair-production. We find 
that 

va=5 X 10-26 cm? sec—! 


about 1/2000 of the rate for the radiative process. Pair-production by 
conversion of the M1 radiation has a similar rate (cf. Rose 1949). 
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ABSTRACT 


The moderation of pile neutrons to very low energies has been investi- 
gated. Liquid hydrogen, liquid hydrogen deuteride and methane at 20°K 
were found to be effective moderators. After an examination of the 
safety aspects, a volume of liquid hydrogen was maintained at the 
centre of a reactor, and the velocity spectrum of neutrons scattered from 
it was measured. Appreciable moderation of neutrons to very low energies 
was observed. By this method it is practicable to increase considerably 
the flux of cold neutrons in a beam. Thus, in the present experiment the 
flux of 10 A wavelength was increased by a factor of 25. 


§ 1. INTRODUCTION 


From the time that nuclear reactors first provided strong sources of 
thermal neutrons, the use of neutrons as a tool for physical research has 
become steadily more widespread. Prominent among the methods 
employed is the scattering of cold neutrons, that is neutrons with 
temperatures less than about 60°K (i.e. with wavelengths >4 A). These 
cold neutrons have two properties of particular value, their low energy 
and their long wavelength. 

Their low energy makes it practicable to measure the energy transfer 
which occurs when neutrons are scattered. This is because the energy 
transfer may be comparable with the initial energy of the cold neutron. 
From a study of energy transfers new information on the thermal 
behaviour of solids and liquids can be obtained. 

Their long wavelength makes cold neutrons suitable for the study of 
long range effects, such as fluctuation phenomena, which are of particular 
interest near the critical point of a phase transition. In addition neutrons 
of sufficiently long wavelengths do not suffer Bragg reflection. Therefore 
they are often used in experiments where Bragg reflection would yield a 
large background of unwanted scattering. 

A third factor, which is likely to assume great importance as the field 
develops, is that a continuous spectrum of cold neutron wavelengths is 
available. Thus, since cold neutrons are obtained by passing a neutron 


+ Communicated by the Authors. 
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beam from a reactor through some form of velocity selector (see e.g. 
Hughes 1953), it is a simple matter to perform an experiment with 
neutrons of various wavelengths. 

As the velocity spectrum of neutrons from a reactor is very nearly a 
Maxwell distribution corresponding to the pile temperature, the flux of 
cold neutrons is relatively small; for example only about 1% of the 
neutrons in a beam from the Harwell reactor BEPO have wavelengths 
greater than 4.4. The fact that only low fluxes are available has limited 
the use of cold neutrons. If, however, one could bring the neutrons into 
thermal equilibrium with a moderator at low temperature, thus shifting 
the whole Maxwell distribution to lower temperatures, the cold neutron 
flux would be greatly increased. 

The use of cooled moderators with cyclotrons is well known (e.g. Sutton 
et al. 1947, and Squires and Stewart 1955); in particular Sutton et al. 
investigated the use of liquid hydrogen as a moderating material. In this 
application, however, it is possible to cool a large volume of moderator, 
whereas in a reactor the volume which could be used would be severely 
limited by the need to situate the moderator in a high neutron flux, Le. 
near the core of the reactor, where the space available for apparatus is 
small. In most materials there would be little moderation to very low 
energies in the volume available since cross sections for inelastic neutron 
scattering fall as the neutron energy becomes comparable with the atomic 
binding energy of the moderator. If a material is to act as a low 
temperature moderator, it must have a high free-atom cross section for 
neutron scattering, a weakly bound structure and also a small absorption 
cross section. 

In a preliminary experiment (§ 2) we showed the feasability of moder- 
ating neutrons to low energies. Of the moderators found to be effective 
liquid hydrogen was the most convenient to use. Therefore, after 
studying the safety of putting liquid hydrogen in a reactor (§ 3), we 
maintained a volume of liquid hydrogen at the centre of the BEPO reactor. 
The velocity spectrum of neutrons scattered from the hydrogen was 
measured. This confirmed that liquid hydrogen was effective in producing 
large fluxes of cold neutrons. 


§ 2. PRELIMINARY EXPERIMENTS 


In the absence of data on cold neutron cross sections, the material which 
will most efficiently moderate neutrons to low energies and the form in 
which it should be used are not obvious. A simple experiment was 
therefore carried out to compare the moderating efficiencies of different 
materials. 

2.1. Apparatus 

The cryostat used for these experiments was similar to the apparatus 
later used in the reactor (§ 4) and is shown in fig. 1. The moderator 
under investigation was contained in a chamber 7-5 cm in diameter and 
10 cmlong. This was cooled by contact with a second chamber containing 
a coolant. The chambers were enclosed in a vacuum case. 
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The cryostat was surrounded by wax in form of a cylindrical block 
45cm diameter and 60cm long. In the wax, four 250 mc polonium— 
beryllium neutron sources were buried at evenly spaced points about the 
cryostat on an 18 cm diameter circle. This circle was 2:5 cm above the 
level of the bottom of the moderator chamber. Thus neutrons reaching 
the moderator chamber had traversed at least 4m of wax. Measure- 
ments with a BF, counter at the position of the moderator chamber 
confirmed that the fluxes of thermal and epi-thermal neutrons falling on 
the chamber were in roughly the same ratio as in a reactor. 


Fig. 1 
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Apparatus used in preliminary experiments. 


Neutrons scattered by the moderator through a collimator formed an 
external beam. Collimation was such that only neutrons coming from 
the direction of the moderator reached the neutron detector placed below 
the wax block. 

Because the beam intensity was very low we could not use a neutron 
velocity selector, and therefore we could only measure an average neutron 
energy for the beam. This was deduced from the attenuation of the 
beam by a gold sheet, } mm thick. On the assumption that the distri- 
bution of neutrons was Maxwellian, the temperature of this distribution 
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was calculated ; we call this the ‘ apparent beam temperature’. If the 
velocity distribution is not Maxwellian (and later measurements, § 5, 
show that in general it is not) this temperature has no absolute significance, 
but it does serve to indicate average neutron velocity and can be used to 
compare the efficiency of different moderators. 

With the moderator chamber filled with pentane at 293°K we verified 
that the apparent beam temperature was correct, viz. 294-+10°K. 


2.2. Results 


Extensive measurements were made on liquid hydrogen since this 
appeared to be the most promising moderator. 

First, the variation of beam intensity with the depth of hydrogen in the 
cryostat was measured. It was found that the intensity was almost 
independent of the depth of hydrogen providing this exceeded 2 cm. 

Figure 2 shows the apparent temperature of a beam scattered from 
liquid hydrogen at 20°K as a function of the relative proportions of ortho- 
and para-hydrogen. The moderation efficiency is seen to vary little over 
the range: 10% to 75% of ortho-hydrogen. 


Table 1. Apparent Beam Temperatures obtained with 
Various Moderators (in °K) 


Temperature of 
moderator 90 77 20 4 


°. 


Moderator 

Pentane — 135+ 14 121215 — 
Methane T3525 (Bybse  7/ (56:5 6548 
Methyl alcohol — (hiss 1G) Vip ayse 17 — 
Hydrogen deuteride — — 77+ 6 = 


Figure 3 shows the effect of varying temperature with hydrogen of fixed 
composition, viz. 45% ortho-hydrogen. The point at 4°K refers to 
solid hydrogen cooled by liquid helium ; the other points refer to liquid 
hydrogen boiling at different pressures. 

In table 1 are given results obtained with other moderators. It is seen 
that only hydrogen deuteride and methane are comparable with hydrogen 
in their ability to moderate neutrons to very low energies. The efficient 
moderation in methane is presumably due to the excitation of rotational 
levels, for methyl alcohol, which is similar in structure but has a greater 
molecular moment of inertia, is much less efficient. 


2.3. Conclusion 
Though liquid hydrogen at various temperatures, hydrogen deuteride 
and methane give similar average beam temperatures, the actual velocity 
spectra may differ considerably in detail ; in particular their yields of cold 
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neutrons may differ. Of the moderators found to be efficient however, 
liquid hydrogen at about its normal boiling point is the easiest to use, and 
so apparatus was designed which would permit a sample of liquid 
hydrogen to be maintained at the centre of BEPOs. 
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Apparent temperature of a neutron beam scattered from liquid hydrogen at 20°K. 
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TEMPERATURE OF HYDROGEN IN °K 


Apparent temperature of a neutron beam scattered from hydrogen of fixed 
composition at varying temperatures. 


+ It has been pointed out by J. Wright (private communication) that methane 
would not be a convenient moderator to use in a reactor, since under irradia- 
tion it undergoes polymerization and evolves free hydrogen, 
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§ 3. SaretTy ASPECTS} 


The presence of some 200 cm* of liquid hydrogen inside a reactor 
constitutes a potential explosion hazard. Besides taking all the normal 
precautions in using liquid hydrogen, we considered it necessary to show 
that even in the event of an explosion there would be ‘no damage to the 
reactor. 

The low temperature apparatus inside the pile was to be enclosed in a 
vacuum case of magnesium—zirconium alloy with walls 6mm thick. To 
test that this would contain any explosion occurring in the apparatus, 
explosions were initiated in an identical tube buried vertically in the 
ground. With the tube filled with the optimum explosive mixture of 
gaseous hydrogen and air, explosions produced shock waves with peak 
pressures of 27 atmospheres at the closed end of the tube. 

Liquid hydrogen explosions were then investigated. It proved very 
difficult to ignite mixtures of liquid hydrogen and solid air. We therefore 
placed separate Dewar vessels of liquid hydrogen and liquid air at the 
bottom of the tube, filled the tube with hydrogen gas and air and 
initiated explosion in the gas. Even when we exploded 300 cm? of liquid 
hydrogen with the optimum amount of liquid air the tube was only 
slightly distorted at the point of explosion ; the peak pressure at the end 
of the tube was 77 atmospheres. The actual vacuum case used in the 
pile was tested repeatedly to a pressure of 47 atmospheres before assembly. 

It is difficult to see how an explosion of the 200 cm? of liquid hydrogen 
inside the pile could occur during an experiment with anything 
approaching the optimum amount of air. A sudden large leakage of air 
into the vacuum space would cause most of the hydrogen to boil off 
before much air had condensed. To prevent a slow leakage causing an 
accumulation of solid air, the vacuum case was not pumped continuously, 
so that the helium and neon content of the air would cause a noticeable 
deterioration in the vacuum before much solid air had collected. 


§ 4. Bepo Liguip HypRoGEN MODERATOR 


An apparatus was installed in BEPo to maintain a volume of liquid 
hydrogen near the centre of the reactor. The chamber containing the 
hydrogen was at the centre of a hole traversing the reactor, and was fed 
with liquid hydrogen from the north face where all the low temperature 
equipment was situated. Neutrons scattered from the hydrogen and 
passing along the hole to the south face of the reactor formed a beam 
available for experimental work. 


4.3. Low Temperature Equipment 
The general lay out of the low temperature equipment is shown in fig. 4, 
and details of apparatus inside the pile in fig. 5. The liquid hydrogen 
moderator was contained in an aluminium can K, 7-5 cm diam. and 


+ Our thanks are due to Mr. J. E. Uppard and Dr, W. G. P. Lamb of A.W.R.E, 
Foulness for their help with these tests, 
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6-25 cm long, inside which was a cooler L with aluminium fins. The 
aluminium tubes (5mm in diam.) supplying the liquid hydrogen were 
supported by wires from a set of cages contained inside the vacuum case 
of magnesium—zirconium alloy (§ 3). This case was 6 m long and extended 
to the pile face. 


Fig. 4 


Lay-out of apparatus used to maintain liquid hydrogen in the BEPo reactor: 
A, Hydrogen gas regulator and relief valve; B, 68 litre liquid hydrogen 
container ; C, moderator chamber; D, control valves; E, duct to 
atmosphere ; F, condensing cylinder and charcoal cleaner; G, relief 
valve (4-7 atmospheres) ; H, oil diffusion pump ; I, ionization gauge : 
J, carbon resistance thermometers. 


Fig. 5 


Details of apparatus within reactor core: K, aluminium can; _ L, cooler ; 
M, beryllium block used in some experiments ; N, evacuated flight tube. 


Neutrons escaping from the reactor down to the vacuum case would 
consitute a radiation hazard. Blocks of neutron absorbing material were 
therefore placed in the vacuum case near the reactor face. 

Two separate hydrogen systems were used. The moderator chamber 
was kept near 20°K by a flow of liquid hydrogen into the cooler L from a 
68 litre storage Dewar B. A constant pressure was maintained over the 
liquid hydrogen in B, and the flow through the circuit controlled by two 
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valves D at the outlet. After the initial cooling the lowest flow we were 
able to use was 2 litres of liquid hydrogen an hour, though mechanical 
failure of the supports of the pipelines and the chamber, later resulted in 
an increased hydrogen consumption. However the system could be used 
satisfactorily even when its consumption was as high as 6 litres an hour. 
The cold gas emerging from the cooler was used to cool a radiation shield 
surrounding the liquid hydrogen line, thus reducing the heat inflow to it. 

The main volume of liquid hydrogen K was condensed into the moderator 
chamber from a cylinder through a charcoal cleaner, and kept at a pressure 
of 5 atmospheres. At this pressure the boiling point is 27°K ; by keeping 
the chamber below this temperature boiling inside the chamber, with 
consequent fluctuations in the cold neutron flux, was prevented. The 
low temperature apparatus could be left unattended in this state for many 
hours. When the liquid hydrogen supply ran out, the hydrogen in the 
moderator chamber simply escaped from the relief valve G. In some of 
the experiments a beryllium block M was attached to the face of the 
moderator chamber (see § 6). 


4.2. Measurement of the Neutron Spectrum 


Neutrons scattered out of the south face of the moderator chamber 
could leave the pile along an evacuated flight tube (N in fig. 4), 9 cm in 
diameter, thus forming an external beam. The flight tube extended to 
the counting equipment, 6 m from the south face of the reactor (i.e. 12 m 
from the moderator). Collimators of boric acid and paraffin wax reduced 
the beam diameter to about 4cm and ensured that the beam consisted 
only of neutrons coming from the direction of the liquid hydrogen. 

The velocity spectrum of the neutrons in the beam was measured with 
the Harwell slow chopper (Hgelstaff 1954), using a BF, neutron detector. 
As the chopper is constructed of cadmium it does not prevent neutrons 
with wavelengths shorter than 0-34 reaching the detector. These 
neutrons would cause a high background count. By passing the beam, 
through a block of beryllium, 23cm long, cooled in liquid nitrogen 
virtually all neutrons with wavelengths shorter than 4 4 were removed 
without serious attenuation of the neutrons of longer wavelength. 


§ 5. Resuuts 


Figure 6 shows the velocity spectrum of the neutron beam obtained 
when the moderator chamber was filled with liquid hydrogen (about 50% 
ortho-hydrogen) at 25°K. The rapid fall in counting rate at a neutron 
time-of-flight of ~1000 y.s/m is caused by the beryllium filter. The error 
bars indicate the statistical uncertainty of the results. Also plotted in 
fig. 6, on the same scale, is the spectrum obtained with the low temperature 
moderator replaced by graphite blocks at the pile temperature (130°c), 
all other conditions remaining unchanged, This second spectrum there- 
fore represents the sort of beam normally used as a source of cold neutrons, 
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In plotting fig. 6 no allowance has been made for the fact that the 
detector efficiency, the transmission of the chopper and the transmission 
of the beryllium all depend on neutron velocity ; the spectra therefore 
differ somewhat from the true velocity distributions. At any given 
velocity, however, the ratio of the neutron flux obtained with hydrogen 
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to that obained with graphite is independent of these effects. This ratio 
has been plotted in fig. 7. The curve (which is now independent of the 
counting system) shows the increase in cold neutrons flux obtained using 
liquid hydrogen instread of a moderator at pile temperature. 

By removing the beryllium filter we extended the measurements of the 
neutron spectrum to higher velocities. The results are shown in fig. 8. 
The need to correct for a large background of fast neutrons now makes the 
statistical accuracy of the results inferior to that of fig. 6, and increases 
the possibility of systematic error. However, it can be seen from the 
shape of the curve that, using this amount of liquid hydrogen, moderation 
to hydrogen temperatures is incomplete. 


Fig. 8 
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Spectrum of neutrons scattered from liquid hydrogen. (No beryllium filter in 
the beam.) 


It should be noted that in measuring this spectrum the angular velocity 
of the slow chopper was twice that used in measuring those shown in 
fig. 6. The consequent change in the effective transmission of the chopper 
is the cause of the slight difference in appearance, over the same energy 
range, of the two curves obtained for neutrons scattered from hydrogen. 

Since a larger sample of moderator could not be used an attempt was 
made to increase the efficiency of the moderator by attaching a disc of 
beryllium, 4 cm thick (M in fig. 5), to the face of the moderator chamber. 
It was hoped that this would reflect fast neutrons leaving the hydrogen, 
thus keeping their flux high in the moderator, whilst the escape of cold 
neutrons would not be impeded. However, the spectrum obtained 
differed little from that shown in fig. 6. 


§ 6. Discussion 


In the absence of data on cold neutron cross sections for liquid hydrogen 
one cannot calculate what neutron velocity spectrum to expect. If 
moderation were complete in the liquid hydrogen the velocities would 
have a Maxwell distribution corresponding to 25°K ; this would have a 
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peak at 1100 us/m on a time-of-flight scale. It is clear from our results 
that only partial moderation has been achieved. The beam consists of 
a broad slowing down spectrum with a low temperature component 
approximately in thermal equilibrium with the moderator, though even 
this component is not exactly Maxwellian in form (fig. 8). 

The first column in table 2 shows the observed gain in neutron flux 
using liquid hydrogen as moderator instead of graphite at pile temperature. 
The second column shows the gain which would be obtained if moderation 
were complete. 


Table 2. Gain in Neutron Flux using Liquid Hydrogen in place of 
Graphite at 130°c 


Neutron wavelength | Observed gain | Gain if moderation complete 


504A 60 


THA 135 
10-04 2! 180 


As a small proportion of the neutrons attain thermal equilibrium with 
the moderator, it is possible that cooling the hydrogen with liquid helium 
might bring some neutrons into thermal equilibrium at 4°x, thus further 
enriching the flux of very long wavelength neutrons, without appreciably 
altering the apparent temperature of the beam as a whole (see § 4). 


§ 7. CONCLUSION 


It has been shown that it is possible and safe to use a low temperature 
moderator, placed in the core of a reactor, to enrich the flux of cold 
neutrons in a neutron beam. The large factors of gain which have been 
obtained with only 200 cm? of moderator make the method very attractive. 
Since the theoretical limit is much higher than the gain observed, the 
method warrants further development. 


ACKNOWLEDGMENTS 


The authors wish to express their thanks to Mr. T. L. Schofield and 
Mr. D. H. C. Harris for their help in assembling and running the apparatus, 
and to Messrs. B. A. Booty, D. H. C. Harris and L. V. Lewis for assisting 
with the measurements. 

REFERENCES 

EaexstaFrr, P. A., 1954, J. nuclear Energy, 1, 57. 

Huauss, D. J., 1953, Pile Neutron Research (Cambridge, Mass.: Addison- 
Wesley, Publishing Co.). 

Seuss, G. L., and Srewart, A. T., 1955, Proc. roy. Soc. A, 230, 19. 

Sutton, R. B., Hatu, T., ANDERSON, E. E., Briper, H. §., DeWrRe, he Wa 
Lavareur, L. S., Lone, E. A., Snyper, T., and Witi1aMs, R. W., 
1947, Phys. Rev., 72, 1147. 


3R2 


[ 928 ] 


The Thermal and Magnetic Properties of Ytterbium Ethyl 
Sulphate between 20°K and 1°K Tf 


By A. H. Cooxz, F. R. McKim, H. MEYER and W. P. WoLF 
The Clarendon Laboratory, Oxford University 


[Received May 13, 1957] 


ABSTRACT 

The susceptibility of a single crystal of ytterbium ethyl sulphate has 
been measured between 20°K and 1°K. Below 4°K the results obeyed the 
relations y,=1-08/7' per mol, and y,=0-061, per mol. These measure- 
ments served to identify the ground doublet state of the ytterbium ion 
in this case as |J=$, J,=+3). The magnetic specific heat was measured 
at five temperatures between 3°K and 1°K by the relaxation method. The 
results were given by OT?/R=4-34 x 10-4, in agreement with calculations 
based on dipole-dipole and hyperfine interactions alone. The spin 
lattice relaxation time + was also measured. At 3°K, 75x 1074 sec and 
was very temperature dependent. 

Susceptibility measurements above 4°K indicated the population of 
excited states, which could not be unambiguously identified. Their 
separation from the ground state was at least 57°K. 


§ 1. INTRODUCTION 


THE rare earth ethyl sulphates form an isomorphous series of salts of the 
general formula M(C,H,SO,); 9H,O where JM is the trivalent rare earth 
ion. They crystallize readily throughout the series in crystals of hexa- 
gonal symmetry. The crystal symmetry results in magnetic symmetry, 
the magnetic susceptibility at any temperature being specified completely 
by principal values parallel and perpendicular to the hexagonal axis of the 
crystal. The crystal structure determination made by Ketelaar (1937) 
shows that all the rare earth ions are magnetically equivalent and the bulk 
properties are therefore closely related to the properties of the individual 
ions. The relation is particularly simple since the crystals are magneti- 
cally so dilute that interaction effects between the magnetic ions are 
comparatively unimportant at temperatures above 1°K. Because of 
these properties, the ethyl sulphates serve as a convenient model series 
for investigation of the magnetic behaviour of rare earth salts. 

Magnetic susceptibility measurements have been made on a number of 
members of the series (de Haas et al. 1933, Fereday and Wiersma 1935) 
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and several have been investigated by measurement of the magnetic 
rotation of the plane of polarized light (Faraday effect) (Becquerel et al. 
1936, 1937, 1938, van den Handel 1940). The results of paramagnetic 
resonance investigations have been reviewed by Bleaney and Stevens 
(1953) and by Bowers and Owen (1955). However, prior to the measure- 
ments described in this paper, no susceptibility measurements had been 
made on the last member of the series, ytterbium ethyl sulphate. Bleaney 
and Scovil (1951) had reported that no paramagnetic resonance could be 
obtained at 20°K and Cooke and Park (unpublished) found that there was 
still no resonance down to 4°K. 

The magnetic properties of the salts depend on the electric crystalline 
field in which the paramagnetic ion is situated. Following a suggestion 
of Ketelaar, Elliott and Stevens (1952, 1953 b) proposed a field of pre- 
dominantly Cy, symmetry, which could be completely specified in any 
particular case in terms of four parameters. A field of this symmetry has 
two effects on any manifold J representing a degenerate state of a free 
ion with an odd number of electrons. The manifold will be split into a 
series of doublets having eigenstates |--/,) where z is the crystal symmetry 
axis, and there will also be admixtures between states with J, differing by 
6 within the manifold. The ground state 4f!° ?F.). of the free Yb** ion 
will therefore be split by a C,,, crystalline field into four doublets. Two 
doublets will be characterized by |-+-4) and |+3), these being pure states ; 
the other two will be admixtures between |+3) and |+3). When a 
magnetic field is applied to the crystal these four doublet levels will be 
split, and, since the condition for there to be an allowed transition between 
two levels is 4./,=0 or +1, the only doublet between whose levels there 
are no allowed transitions is the |3) state. If this doublet were lowest 
then the absence of a paramagnetic resonance spectrum would be 
explained. This was first pointed out by Elliott and Stevens, who found 
that a wide variety of crystal field parameters led to a scheme of energy 
levels having this state lowest. In general the crystal field will also give 
rise to admixtures between the lowest manifold and excited states, as in 
the case of cerium (Elliott and Stevens 1952). In most cases these are 
small and for ytterbium they may be neglected when considering the 
ground multiplet, as the next J level is some 10 300 cm~? higher. 

From measurements on cerium ethyl sulphate, Elliott and Stevens 
deduced the numerical values of the crystalline field parameters for this 
salt. They then suggested how the parameters might change from 
member to member along the rare earth series, and showed that the 
available data could be satisfactorily explained, although the extra- 
polation became increasingly uncertain towards the end of the series. 
It would be particularly interesting, therefore, to determine the values of 
the parameters for ytterbium, the last magnetic member of the series, as 
the values for intervening members could then be obtained by inter- 
polation rather than extrapolation. The experiments described here 
served to identify the ground state unambiguously as |+-3) and showed 
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that the next doublet was at least 57°K higher, but for reasons discussed 
below more general conclusions regarding the crystal field parameters 
could not be drawn. 


§ 2, Tot MEASUREMENTS 


The susceptibility measurements were made on a single crystal using 
the mutual inductance apparatus described elsewhere (McKim and Wolf 
1957). The crystal was prepared from ytterbium oxide supplied by 
Johnson Matthey whose spectroscopic analysis indicated a purity of at 
least 99°, the main impurity being samarium oxide (0-5%). Both the 
temperature dependent and temperature independent contributions to 
the susceptibility were measured. At temperatures in the liquid helium 
range the susceptibility in the direction parallel to the crystal axis was 
found to obey the law 


. “05 
No temperature independent susceptibility was detected, the upper limit 
being 0-01 per mol. From measurements perpendicular to the crystal axis 
it was found that 


x1 = 0-061;-- 0-001 per mol 


independent of temperature throughout the helium range. 

Measurements of y,, and y, were also made in the liquid hydrogen range. 
It was found that perpendicular to the crystal axis the susceptibility 
deviated from the constant value obtained at lower temperatures (see 
figure), and that parallel to the axis Curie’s law was approximately 
obeyed but with a Curie constant somewhat less than in the helium region. 
These measurements indicated that y, should be equal to y, at a tem- 
perature of 18-0°K and this made it possible to check the consistency of 
our results by the following very simple experiment. A crystal was 
suspended with its axis horizontal by a fine thread inside a cryostat whose 
temperature could be varied by pumping the refrigerant (liquid hydrogen). 
A large horizontal magnetic field was applied and the crystal was made 
to oscillate about its axis of suspension. The period of oscillation, t, is 
then proportional to (|x — x,|)~1* and is therefore infinite when yy= ,. 
é was measured at temperatures near 18°K, and from a graph of ¢-? against 
T it was found that y,— x, at 18-4°K in good agreement with the value 
expected from the separate measurements of y, and y,. 

Magnetic specific heat measurements were made in the liquid helium 
region by the relaxation method, using the apparatus described by 
Benzie and Cooke (1950). Measurements were made at five temperatures 
between 1° and 3°xK. The specific heat was found to be inversely pro- 
portional to the square of the absolute temperature. In the course of 
these experiments the spin-lattice relaxation time, 7, was estimated at 
each temperature and it was found that above 3°x, 7 becomes so short 
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that valid results could no longer be obtained with our apparatus. Values 
for the magnetic specific heat and also for 7 in the range 1° to 3°K are 


given in the table. 


x1 10? per mol 


Ytterbium ethyl sulphate. 


perpendicular to the hexagonal axis of the crystal. 


Pid bs ce 


Variation with temperature of the susceptibility 


The symbol + 


indicates the point where yy, as determined by the oscillation 
The dotted line is the value of x, below 4-2°K, 


experiment (see text). 


1-09 — 


1-40 — 2-20 x 10-4 
2°16 2510 | 9-37 x 10-° 
2-68 2 LO OT X10 
2-93 GPL ie O03 x 107" 


§ 3. Discussion 


The results of the susceptibility measurements in the helium region 
may be compared with the general form for the susceptibility of an 


C/R per mol 


3-64 x 10-* 


CT*/R per 
mol 


4-32 x 10-4 

4-31 x 10-4 
4-37 x 10-4 
4-36 x 10-4 
4:32 x 10-4 
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assembly of non-interacting ions each in a doubly degenerate ground state, 
ING GP" (1) 
where N is Avogadro’s number, 8 the Bohr magneton, k Boltzmann’s 
constant, g, the spectroscopic splitting factor for the two levels of the 
ground state when the measuring field is in the ¢ direction, and «, is a 
temperature independent term depending on matrix elements to excited 
states. Comparison of this formula with the results obtained gave 
gu=3-40-0-07 %= 0-40-01 
g,=0-00-40-05 a, =0-061- 0-001 

Of the four possible doublet ground states only |+3) has gy,=0. The 
value of g, for this state is equal to 3g;, where g,=8/7, the Landé g factor 
for the 2Fj,. manifold. This gives g,—3-43 in excellent agreement with 
the experimental value. From this we conclude that the ground state 
is indeed |J,=-+3), as suggested by Elliott and Stevens. ; 

The magnetic specific heat measured at the lowest temperatures, when 
all the ions are in the ground state, contains two contributions, one due 
to interactions between neighbouring paramagnetic ions in the crystal 
lattice, the other due to the splitting of the ground state of each ion by 
nuclear effects. If the interaction specific heat be attributed solely to 
magnetic dipole-dipole interaction, which is likely as exchange effects 
are known to be small in the isomorphous neodymium ethyl sulphate 
(Roberts et al. 1953), then by standard methods (Van Vleck 1937, Daniels 
1953) it can be shown that 

Caa__ _gu'B* t zis zi" 

taser Bap 8B +9 R745 | +o(s) - 
where 7,, is the distance between the ions 7 and j and z,; is the component 
of this along the z-axis. In the temperature range of our measurements 
only the first terms in this expansion need be considered, and using the 
values of the lattice sums evaluated by Daniels, eqn (2) reduces to 

2 4p4 
os = ae (532-2)a-6 

when a is the unit cell dimension perpendicular to the hexagonal axis. 
Taking the value af a as 13-814, obtained by extrapolation of this 
quantity for earlier members of the rare earth series (Ketelaar 1937), we 


find 


2 
The other contribution to the specific heat, due to hyperfine interaction 
had to be evaluated in an indirect manner in view of the absence of para- 
magnetic resonance data for this salt. The general spin Hamiltonian 
describing the interaction of an ion in a doubly degenerate electronic 
ground state with both a magnetic field and its nucleus is given by 


H#=PH.g.S+S.A.I ..... . (3) 
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(Abragam and Pryce 1951), neglecting small quadrupole effects. If the 
electric field which has given rise to the doublet has axial symmetry this 
reduces to 


HA =GyPHS,+9,8(AS,+H,S,)+A,S1,+ 41 (Sol gtS,1,)- 
If admixtures from excited states having a J differing from that of the 
ground manifold can be neglected, we have a simple relation between the 
principal values of the gy and A tensors (Elliott and Stevens 1953 a) 
Ay _ Ay 
Mi = 9h 
where y is a constant for the ion independent of the particular crystal 


field in the salt. Thus, for the state |+3), for which g, = 0, A ,=0, and 
the spin Hamiltonian reduces to the simple form 


H = =9PH, S ert (vawSL,. 


The constant y was calculated from the resonance data on ytterbium 
acetate (Cooke and Park 1956) for which transitions within the lowest 
doublet are not forbidden as in the ethyl sulphate. 

In zero field a spin Hamiltonian of the form (3) leads to a specific heat 
contribution 


aw 


Cy The? i 
R =5 (zr) {(yau)?S(S+ DLL +1)} 


(Bleaney 1950). Substituting the value of y obtained from the acetate 
and allowing for the abundances of the two Yb isotopes with non-zero 
spin, we find 
yp 
oe = 3-22 x 10-4. 

The sum of the two contributions yields a total specific heat given by 
CT?/R=4-46 x 10-4, compared with the measured value of 4:34 10-4. 
The close agreement between these confirms that there is very little 
contribution from exchange interaction to the specific heat, which can 
be explained in terms of magnetic dipole-dipole and hyperfine inter- 
actions alone. 

The effect on the low temperature susceptibility of the dipole-dipole 
and nuclear interactions can similarly be calculated in terms of the lattice 
parameters, g value and the constant y. Expressing the correction in 
terms of a power series in 1/7 it can be shown (Daniels 1953) that dipole— 
dipole interaction contributes a term —(0-013/7')x9, where yo is the first 
order term in 1/7’ given by eqn (1) while the effect of nuclear interaction 
enters only in the next approximation (Bleaney 1950). At temperatures 
down to 1°K these corrections are therefore negligible and Curie’s law is 
obeyed. Below 1°k they become increasingly more important, but even 
down to 0:1°K the deviations from Curie’s law will be small. 

The deviation of the susceptibility in the hydrogen region from its 
behaviour in the helium region was interpreted as due to some ytterbium 
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ions being excited into states other than the ground doublet. Population 
of these excited states will also give rise to a Schottky anomaly in the 
specific heat, which at the low temperature end has the form 
OT?2|R=A* exp (ap 

where 4 is the energy of the first excited state, expressed in °K. This 
result is quite independent of the magnetic strength of the excited states 
since the specific heat depends only on the energy difference between the 
states, and not on their quantum. description. Now at the highest tem- 
perature of measurement, 2-93°K, no onset of an anomaly was observed. 
Assuming that a 5°% variation in CT?/R would have been detected, a lower 
limit could be put on 4, namely 4=57°K. Other than this, it was not 
possible to draw any definite conclusions about the identification and 
positions of the excited states. Calculations showed that the suscepti- 
bility results could not be explained with the |+4) state as the first 
excited state and the other two states at much higher energies, whereas a 
reasonable fit was obtained with one doublet, principally |+$) with a 
small admixture of |--3) at about 70°K above the ground state. But the 
possibility of there being two doublets at about the same energy separation 
from the ground state could not be excluded. It would be very desirable 
to have specific heat measurements at temperature above 3°K. If the 
onset of the Schottky anomaly could be observed then it would be possible 
to decide whether one or two excited states were being populated, and to 
obtain some idea of the splitting between these levels and the ground state. 

The susceptibility measurements have shown that at temperatures 
below 4°x ytterbium ethyl sulphate is completely anisotropic. Ifa single 
crystal is demagnetized adiabatically then the subsequent application of a 
large magnetic field perpendicular to the crystal axis will produce no 
heating. It is therefore a better substance than is cerium ethyl sulphate 
for the type of magnetic experiments below 1°K suggested previously by 
Bogle et al. (1951). Measurements of the thermodynamic properties 
below 1°K have recently been made and will be published in due course, 
together with experiments in nuclear alignment. 
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A ctoup chamber evidence for cascade decay of a heavy charged 
unstable particle is reported here. The picture shows two V-events, 
one due to a charged particle and the other due to a neutral particle. 
From a consideration of the coplanarity of the apex of the charged 
V-event and the tracks of the two secondaries of the other and also of 
the angles of emission and minimum ionization of all the decay products 
it appears that a heavy particle of mass not below 1460 m, has undergone 
a cascade decay process according to the mode 


K+ + L++6 +> 7t+n-+212 mev. 


The photograph reproduced in fig. | (Plate), shows that a charged 
particle PO, which is produced in an interaction at P (Cu-plate), traverses 
about 5 cm with 2-3 times minimum ionization and then decays in the 
gas at O. The charged decay product OA starting from O disappears 
towards the rear of the chamber at minimum ionization after passing 
through 29 g of Cu equivalent. At O’ a V°-particle is found to decay 
into two minimum ionizing particles O'B and O'C, the former disappearing 
through the left side-wall and the latter through the front plate of the 
cloud chamber after penetrating 15g of Cu equivalent without 
multiplication. 

Since both the events take place in the gas volume, we exclude the 
possibility of a capture process at O or a pair creation of electrons at O’ 
and proceed to correlate the two events. A stereo-projection shows that 
the decay point O of the charged V-track POA lies in the plane containing 
the two decay tracks O’B and O’C of the neutral V-event within +3°. 
This strongly suggests that the particle PO has a neutral decay product 
OO’ besides the charged one OA, and this neutral particle most probably 
produces the decay event at O’. The angles made by the emitted particles 
with their respective parents in space are shown schematically in fig. 2. 


INTERPRETATION 
Since most of the neutral V-events are either A° or 6° decays, we 
naturally try to identify the V°-event at O’ with one of them. If we 
assume it to be a case of A° decay, it follows that the particle O’B making 
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the smaller angle of 27° with the line of flight of the neutral particle 
should be a proton and the other particle O’C making the angle 65° a 
m-meson. Since O’B is a proton travelling with minimum ionization 
its momentum must exceed 1 Gev/e. Conservation of momenta and 
energies then leads to a minimum mass of 1-7 Gev for the neutral particle 
OO’, which rules out the possibility of identifying it with a 4°. A similar 
calculation with the assumption that it is a case of 6° decaying into two 
pions gives a minimum mass value of 468 Mey for the V°-particle. Here, 
we take the momentum of the vertical secondary O’C which makes the 
larger angle and remains at minimum ionization after passing through 
15 g of Cu, to be at least 180 mev/c at the point of emission. The lower 
limit of 468 Mev for the mass of the neutral particle is not far from the 
known 6°-mass of 493 Mev. We, therefore, conclude that the neutral 
particle OO’ is a 6°-particle. 


Fig. 2 
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The various angles made in space by the different secondaries with their 
respective parents are shown schematically. 


The identification of the charged decay product OA is made in an 
indirect way. Since OO’ makes an angle (20°) which is much smaller 
than that (60°) made by OA with the parent particle PO, we can reasonably 
conclude that the particle OA is considerably lighter than OO’. The 
particle OA is therefore taken to be an L-meson of a minimum momentum 
190 Mev/e at the point of emission corresponding to an almost minimum 
ionization even after passing through 29 g of Cu equivalent. Assuming 
a two-body decay of PO, conservation of transverse momentum requires 
a minimum momentum of 481 mMev/ce for the neutral secondary OO’, 
which has been previously identified with a 6°. Conserving energy and 
longitudinal momentum the minimum mass and momentum which we 
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find out for the primary particle PO are 745 Mev (1460 m,) and 549 Mev/c 
respectively pointing to an ionization of 2-5 times minimum which is 
quite consistent with its observed ionization. 

All cascade decays so far reported (Thompson 1956) have A® as the 
neutral secondary and one leg of the secondary decay event has been 
invariably of more than minimum ionization. In the present picture 
both the secondaries of the secondary decay event at O’ are distinctly 
of minimum ionization and therefore its identification with a A° decay is 
not possible. 

Another possibility is that the 6° decaying at O’ and the K-particle 
decaying at O might have been produced in the same interaction at 2) 
This is very unlikely since the plane containing O’B and O’C makes an 
angle of more than 15° with the point P obtained by joining the three 
visible tracks emitted from that plate. 


LIFETIMES OF THE PARTICLES 


As both the decays have occurred in flight and the decaying particles 
are both produced within the visible region, a lifetime determination is 
possible. The paths traversed by the charged and neutral decaying 
particles are 5-1 and 16+2 cm respectively in the laboratory system 
between production and decay. We, however, know only the lower 
limit of their momenta and this gives an upper limit for the lifetimes as 
2510-1 and 6:5x10-1°sec for the charged and neutral particles 
respectively. From the above considerations it appears that the two 
V-events in fig. 1 (Plate) are actually linked with each other and they 
represent a case of cascade decay of a charged K-particle whose mass is 
above 1460 m, and whose lifetime is below 2-5 x 10-1 second. 
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REVIEWS OF BOOKS 


Theories of Nuclear Moments. By R. J. Burn-Sroyte. (Oxford University 

Press.) [Pp. 88.] Price 8s. 6d. 

In these days of intense specialisation, when papers on any one subject appear 
in many journals and advanced textbooks are invariably expensive, the 
introduction of a new series of inexpensive paper-backed monographs called 
the Oxford Library of the Physical Sciences is most welcome. 

This volume, the first of the series, is mainly concerned with a comparison of 
the predictions of theories based on simple nuclear models with the experimental 
results. It contains an excellent account of moment calculations using the 
different particle and collective models, and the effects of meson currents, 
velocity dependent forces, intermediate coupling and interconfigurational 
mixing are discussed in some detail. A chapter on the experimental methods of 
measuring nuclear spins and moments is included for completeness, and several 
tables of moments, both experimental and theoretical, with notes about the 
relevant models are set out in the text and in an appendix. KES: 


Fatigue in Aircraft Structures. Proceedings of the International Conference 
held at Columbia University, January 30th to February Ist, 1956. Edited by 
A.M. FREUDENTHAL. (New York: Academic Press Inc. ; London: Academic 
Books Ltd.) [Pp. ix + 456.] Price $12.00. 

THis book contains a number of contributions on the effect of fatigue on 

design of aircraft and on cumulative damage under varying stress. In addition 

there are several very interesting review papers, by Wood, Forsyth, Thompson 
and others on the physical processes observed as a fatique crack is developing, 
though it is not made clear whether this work has as yet contributed to the 

development of alloys for practical purposes. N. F. M. 


Progress in Low Temperature Physics. Edited by C. J. Gorrrr. (North 

Holland Publishing Company.) [Pp. xi + 480.] Price 84s. 

THE second volume of ‘“ Progress in Low Temperature Physics ’’, like its 
predecessor, is of very uneven quality and interest, but it possesses enough 
articles of value that one may welcome its appearance. It may however be 
queried whether Low Temperature Physics is nowadays a sufficiently self- 
contained topic to justify a regular examination of its progress. Certainly 
some of the articles here on solid-state matters, such as the one concerned 
with semi-conductors, suffer from too rigid an attention being paid to the 
low temperature aspects, while at the other extreme the article on the rare- 
earth metals hardly seems to qualify for admission into this volume. Of the 
traditional problems of low temperature research, liquid helium is well represen- 
ted by 137 out of the 480 pages, while 13 is all that is allowed to super-conducti- 
vity. Perhaps this bias, already noticeable in the first volume, may be taken 
to reflect the special interests of the editor, whose hand is also understandably 
discernible in the choice of authors. This tendency, normally beneficial or 
at any rate harmless, is rather regrettable when it lends itself to partisanship 
in the current controversy over an improved temperature scale, a matter which 
has already led to acrimony and in which the right is by no means obviously 
on one side or the other. 

Much that is discussed here is of such ephemeral interest, and that only 
to the low temperature specialist, that it is encouraging to find extended 
accounts of topics which in the future may have something to contribute in a 
wider sphere than mere low temperature physics ; the articles by Shoenberg 
on the de Haas-van Alphen effect and by Steenland and Tolhoek on nuclear 
orientation are on this account particularly welcome, and help to redeem the 
yolume of something of the atmosphere of a parish magazine. AL BSP: 


940 Reviews of Books 


The Defect Solid State. By T. J. GRay, D. P. DetwimezrR, D. E. Rass, W. G. 
Lawrence, R. R. West and T. J. Jennines. (New York: Interscience 
Publishers Inc.; London: Interscience Publishers Ltd.) [Pp. vii+511.] 


Tuts book is a gallant attempt on the part of members of the staff of the 
College of Ceramics of Alfred University, New York, to put together the funda- 
mental parts of the physics and chemistry of the solid state which are relevant 
to the branch of technology with which the College is concerned. This includes 
the properties of dielectrics, magnetism, corrosion and in particular reactions 
in solids, Inevitably, through covering so wide a field, the treatment is 
sketchy in places; but none the less this attempt to provide a textbook of 
fundamental theory for a particular industry is of real interest. N. F. M. 


Photographic Sensitivity. Vol. 1. Hakone Symposium. Edited by Shin 
Fujisawa. (Tokyo: Maruzen Co. Ltd.) [Pp. 147.] Price $4.00. 


Unper the title ‘‘ Photographic Sensitivity Vol. 1 ’’ published by the Maruzen 
Company, Tokyo, 1956, are published the Proceedings of a Symposium held 
at Hakone in September 1953. There were some forty attenders, mostly 
Japanese, but including N. F. Mott and F. C. Frank from Bristol and F. Seitz 
from the University of Ilinois. The main subjects covered The Physical 
Properties of Silver Halides and Photographic Emulsions, ranging from lattice 
defects and the orientation of print-out silver with respect to the parent silver- 
halide crystal to the chemical sensitization of emulsions and the mechanism 
of sensitivity of collodion emulsions. This successful venture.was continued 
on two later occasions by The Group of Research Workers on Lattice Defects, 
and it is planned to issue further publications in the same form. Wi bs 5 


Elements of Partial Differential Equations. By I. N. SNEDpon. (London 
and New York: McGraw-Hill.) [Pp. ix+327.] Price $7.50. 


Tuts is an account of methods for finding analytical solutions of partial differ- 
ential equations. There is more to it than would be needed by most under- 
graduates, but the unsophisticated treatment of many different techniques 
should make it extremely useful for theoretical physicists, engineers, and other 
applied mathematicians. There are numerous examples and problems. 


J.M.Z. 


Scientific Uses of Harth Satellites. By James A, vAN ALLEN (ed. by). (Chap- 
man & Hall.) [Pp. x+316.] Price 63s. 


Tuts book contains 33 articles given at a meeting of the U.S. Upper Atmosphere 
Rocket Research Panel in January, 1956. The subjects covered are satellite 
orbits, methods of tracking, instrumentation and measurements using satellites, 
both from the ground and from instruments in the satellite. The subjects 
whose investigation is considered include the figure of the earth, air density 
and composition, meteorology, meteorites, solar radiation, cosmic rays, 
terrestrial magnetism, the ionosphere and the aurora. 

The articles are of considerable interest. They represent practicable experi- 
ments most of which we may expect to see carried out during the next twenty 
years. One of the most ambitious projects is to fly a mass spectrometer in a 
satellite. E. C. B. 


[The Editors do not hold themselves responsible for the views 
expressed by thei correspondents. | 
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Left view 


Right view 


Photograph showing the decays of a charged particle at O and a neutral particle 
at O’. The plarie containing the secondaries O’B and O’C of the V°-event 
at O’ passes through O within +3°. It is interpreted that the mass of 
the particle PO is above 1460 m, and it decays into an L-meson OA 
and a 6° represented by OO’ which subsequently decays into two pions 


O’B and O'C. 


